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This new calibration atlas is based on 
frequency rather than wavelength cali- 
bration techniques for absolute refer- 
ences. Since a limited number of 
absolute frequency measurements is 
possible, additional data from alternate 
methodology are used for difference 
frequency measurements within each 
band investigated by the frequency 
measurements techniques. Data from 
these complementary techniques include 
the best Fourier transform measure- 
ments available. Included in the text re- 
lating to the atlas are a description of 
the heterodyne frequency measurement 
techniques and details of the analysis, 
including the Hamiltonians and least- 
squares-fitting and calculation. Also in- 
cluded are other relevant considerations 
such as intensities and lincshape 
parameters. A 390-entry bibliography 
which contains all data sources used 
and a subsequent section on errors con- 
clude the text portion. 

The primary calibration molecules 
are the linear triatomics, carbonyl sul- 
fide and nitrous oxide, which cover por- 
tions of the infrared spectrum ranging 
from 488 to 3120 cm -1 . Some gaps in 
the coverage afforded by OCS and N2O 
are partially covered by NO, CO, and 
CS2. An additional region from 4000 to 
4400 cm -1 is also included. 



The tabular portion of the atlas is 
too lengthy to include in an archival 
journal. Furthermore, different users 
have different requirements for such an 
atlas. In an effort to satisfy most users, 
we have made two different options 
available. The first is NIST Special 
Publication 821, which has a spectral 
map/facing table format. The spectral 
maps (as well as the facing tables) are 
calculated from molecular constants 
derived for the work. A complete list of 
all of the molecular transitions that 
went into making the maps is too long 
(perhaps by a factor of 4 or 5) to in- 
clude in the facing tables. The second 
option for those not interested in maps 
(or perhaps to supplement Special Pub- 
lication 821) is the complete list (tables- 
only) which is available in computerized 
format as NIST Standard Reference 
Database #39, Wavelength Calibration 
Tables. 
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1. Introduction 

The primary purpose of this work is to provide an 
atlas of molecular spectra and associated tables of 
wavenumbers to be used for the calibration of in- 
frared spectrometers. A secondary purpose is to 
furnish a detailed description of the infrared het- 
erodyne frequency measurement techniques devel- 
oped for this work. Additionally, we provide a 
bibliography of all the measurements used in pro- 
ducing these tables, as well as to provide a descrip- 
tion of how those measurements were combined to 
calculate energy levels, transitions, and uncertain- 
ties. We also provide useful related information 
such as line intensities, pressure broadening coeffi- 
cients, and estimates of pressure shifts of spectral 
lines. This book does not include an exhaustive list 
of all the weaker transition frequencies currently 
available, especially for the less abundant molecular 
species. Such a list, containing over 10 000 transi- 
tions for OCS, is available, however, in computer- 
ized format as NIST Standard Reference Database 
39. To put this work in proper perspective, some 
background, philosophy, and the status of existing 
atlases are discussed in the following sections. 
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Over the last 35 years (since the work of Downie 
et al. 3 [1.1]) several compilations of infrared ab- 
sorption spectra intended for the calibration of in- 
frared spectrometers have appeared. Two compi- 
lations have been published by the International 
Union of Pure and Applied Chemistry (IUPAC) 
[1.2,1.3] and contain sections pertaining to the cali- 
bration of fairly low resolution (0.5 cm -1 ) instru- 
ments. Other compilations [1.4-1.7] and other 
sections of the IUPAC compilations [1.2,1.3] were 
devoted to data intended for the calibration of high 
resolution instruments (resolution better than 0.5 
cm -1 ). This book falls into the latter category. Of 
the earlier compilations, only the work of 
Guelachvili and Rao [1.7] provides calibration data 
that are consistently more accurate than ±0.01 
cm" 1 . 

A number of commercially available infrared 
spectrometers are capable of recording spectra 
with a resolution of 0.06 cm -1 or better. For the 
calibration of such spectrometers, one needs cali- 
bration data with an uncertainty smaller than 1/ 
30th of the spectral linewidth. For a resolution of 
0.06 cm -1 this means an accuracy of 0.002 cm -1 . 
For Doppler-limited resolution at room tempera- 
tures, a calibration accuracy of 0.0002 cm -1 or bet- 
ter would be desirable. The present work responds 
to these needs, even though it was recognized that 
state-of-the-art instrumentation (for sub-Doppler 
measurements, for example) could use even more 
accurate calibration data. 

Most high resolution spectrometers are not ca- 
pable of broad frequency scans. For the most accu- 
rate measurements the calibration should be 
applied to each spectral scan. This requires that at 
least one and preferably several calibration points 
be available within the scanning range of the spec- 
trometer. For tunable diode laser spectrometers, a 
single scan may cover only 0.5 cm -1 , while for a 
Fourier transform spectrometer (FTS) a high reso- 
lution scan may cover hundreds of wavenumbers. 
This means that for many purposes calibration 
standards should be no more than 100 cm -1 apart, 
while for other purposes standards no more than 
0.5 cm -1 apart are required throughout the in- 
frared region of interest. Our work here presents a 
compromise between these two requirements. 

In our opinion, all previous compilations of in- 
frared wavenumber standards have suffered from 
the lack of a consistent effort to draw together a 
number of experimental measurements to arrive at 



3 The references in Sec. 1 are listed at the end of the section as 
[1.1], [1.2], etc. 



a well determined set of molecular energy levels 
which could be used to determine frequency (or 
wavenumber) standards for a number of bands 
throughout the infrared spectral region. The 
present compilation provides a model for produc- 
ing such frequency or wavenumber calibration 
data. 

It is appropriate at this point to define some of 
the terminology used in this book clarifying an im- 
portant distinction: in some places emphasis is 
placed on the difference between frequency mea- 
surements and wavelength measurements. All mea- 
surements that can be reduced to counting 
frequencies or frequency differences are frequency 
measurements, while all measurements that are re- 
ally comparisons of wavelengths (or wavenumbers), 
including counting fringes, are wavelength mea- 
surements. Fourier transform measurements, for 
instance are truly wavelength measurements even 
though counting (of fringes of a laser beam passed 
through the spectrometer, for example) may be in- 
volved. 

Most measurements of infrared absorption spec- 
tra are made using grating or Fourier transform 
spectrometers. These measurements are truly 
wavelength measurements. They rely for calibra- 
tion on measuring the difference in wavelength of 
some calibration feature and the feature of the 
spectrum to be determined. For FTS instruments, 
the position of the moveable mirror must be deter- 
mined, while for grating instruments one is essen- 
tially calibrating the grating angle and the spacing 
of the grooves of the grating. Very often these in- 
struments use more than one beam and compare 
the wavelength of a calibration beam with the 
beam carrying the spectra of interest. With FTS 
instruments a helium-neon laser beam is often 
used to monitor the mirror position. In other cases 
a single beam is used, but one cannot be certain 
that the beam properties are independent of wave- 
length. With all wavelength measuring instruments, 
it is extremely important that the calibration light 
beam follow precisely the same path through the 
spectrometer as the beam of light being measured. 
Problems of wavelength-dependent diffraction ef- 
fects become very important for measurements in- 
tended to approach or exceed one part in 10 7 . One 
of the most insidious problems with wavelength 
measurements is the risk (or ease) of incurring sys- 
tematic errors and the virtual impossibility of de- 
tecting them. 

The great advantage of frequency measurements 
is that they simply depend on frequency counting 
techniques. Over the years considerable attention 
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has been given to electronic techniques of fre- 
quency counting, and simple, accurate, well cali- 
brated devices are readily available. The accuracy 
of frequency techniques does not depend on 
whether or not the beams are perfectly collinear or 
have parallel wavefronts. If a countable signal 
(25 dB signal-to-noise ratio) is obtained, it will be 
correct. The frequency of light does not depend on 
the medium, nor does it depend on the angle from 
which it is viewed. 

Over the past decade we have striven to provide 
infrared heterodyne frequency measurements on 
the transitions and energy levels of several simple 
molecules, which are good candidates for fre- 
quency calibration standards. OCS, CO, and N 2 
were chosen because they are stable, safe to han- 
dle, easily obtainable, and well documented in 
terms of good measurements reported in existing 
literature. Furthermore, they are particularly 
amenable to the accurate calculation of energy lev- 
els from a relatively simple Hamiltonian. 

For these molecules a least-squares fit of many 
transitions (over 3000 OCS lines for example) has 
permitted the determination of all of the lower en- 
ergy separations, using frequency differences re- 
ferred to the primary cesium frequency standard. 
Due to statistical improvements from a large data 
base, transition frequencies between these energy 
levels can be calculated with greater accuracy than 
any single measurement with its attendant random 
errors. 

Although we have had to use some FTS (wave- 
length) measurements to help define certain higher 
order rotational constants, for the most part the 
energy levels were determined horn frequency mea- 
surements because they are less susceptible to un- 
known systematic errors than are wavelength 
measurements. Particular importance was attached 
to estimating the uncertainties in the transition 
wavenumbers; see the discussion in the chapter on 
errors (Sec. 4). 

As with any good calibration standard, a number 
of different measurements were used in determin- 
ing the energy levels and transition frequencies. 
However, only a few laboratories have used fre- 
quency measurement techniques in the infrared re- 
gion and very few of the more accurate 
sub-Doppler frequency measurements have been 
made, so it is somewhat premature to claim the 
level of accuracy that we desire for infrared stan- 
dards. Nevertheless we are encouraged by the con- 
vergence of different FTS measurements on the 
same values for the band centers as frequency mea- 



surements. Our publication of this atlas at this time 
is dictated by the need to provide good calibration 
data now, rather than await the arrival of a perfect 
atlas. 

Of course the combination of OCS, CO, and 
N2O gases is insufficient to provide calibration data 
everywhere within the infrared, so we have had to 
provide heterodyne measurements on other molec- 
ular species, such as CS2 and NO in order to fill 
some of the gaps. The user will still note that many 
gaps remain in the coverage of these tables. To 
some extent these gaps may be filled by using the 
data provided by the compilation of Guelachvili 
and Rao [1.7]. We also expect that future measure- 
ments will provide calibration data where none are 
currently available. 

Since most workers are more comfortable with 
calibration data given in wavenumber units, the ta- 
bles in this book are given in wavenumbers (cm -1 ) 
even though the values were primarily determined 
from frequency measurements. The conversion 
from frequency units to wavenumber units was 
made by using the defined value of the velocity of 
light, c =299 792 458 m/s. Since the tables are given 
in wavenumbers, we often use the terms wavenum- 
ber and frequency interchangeably in the text, but 
the term wavelength is reserved for quantities de- 
termined by wavelength measurements and must 
not be confused with frequency measurements. 
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2. Techniques Used for Infrared Hetero- 
dyne Frequency Measurements 

2.1 Preliminary Considerations 

At the present time the only frequency measure- 
ments that have been made on infrared absorption 
spectra of molecules treated in this book are those 
made in the NIST Boulder Laboratories [5.73], 4 the 
Harry Diamond Laboratories [5.88], the University 
of Lille [5.126], and University of Bonn [5.321]. 
The measurements of the N2O laser transition fre- 
quencies made at the National Research Council 
of Canada [5.183] were also used here. This chap- 
ter describes the techniques that have been used in 
the NIST Boulder frequency measurements in or- 
der to familiarize the reader with the techniques 
that have evolved. We hope that such familiarity 
will give greater confidence in the accuracy of the 
results. 

The first heterodyne frequency measurements 
on carbonyl sulfide (OCS) were made in the NIST 
Boulder Labs. A frequency-stabilized CO2 laser 
served as the local oscillator for these heterodyne 
measurements. For this chapter, the term local os- 
cillator is reserved for the fixed-frequency oscilla- 
tor which carries the reference frequency 
information (near the frequency of the molecular 
transition to be measured) to the mixer or hetero- 
dyne detector. For each OCS measurement, the 
frequency of a tunable diode laser (TDL) was 
locked to the peak of a selected (OCS) absorption 
line with an assigned uncertainty, Sfiock, given by 
± 1/24 ifDopp/SNR, where A n>o PP is the full Doppler 
width at half maximum and SNR is the signal-to- 
noise ratio of the first derivative lock signal. The 
frequency of the locked TDL was compared with 
the frequency of the C0 2 laser frequency standard 
by mixing the TDL and CO2 laser radiation in a 
fast (1 GHz, 3 dB bandwidth) HgCdTe detector. 
The resulting difference-frequency beatnote was 
measured with the aid of a spectrum analyzer and a 
marker oscillator (a conventional oscillator whose 
frequency was tuned to the center of the beatnote). 
The frequency of the OCS transition was the sum 
of the CO2 laser frequency and the appropriately 
signed beatnote frequency. 

Similar sets of measurements have also been 
made in which the TDL was heterodyned against a 
CO laser local oscillator (or transfer oscillator). In 
these cases the frequency of the CO laser was 
simultaneously measured relative to a frequency 

4 Rather than repeat references given in the bibliography 
(Sec. 5), they are referenced as [5.XX]. 



synthesized from combinations of CO2 laser stan- 
dard frequencies [5.94]. In a third type of measure- 
ment, a tunable color center laser (CCL) was the 
local oscillator. In this case, the CCL (which was 
locked to the transition of interest) was hetero- 
dyned directly with a frequency synthesized from 
CO2 laser standards [5.304]. All three techniques 
are discussed in detail later in this chapter. 

The determination of the spectroscopic con- 
stants for a particular band of a selected molecule 
was based on a number of measurements such as 
those described above. Some preliminary consider- 
ations for selection of the particular transitions to 
be measured were as follows. A review of the set of 
transitions of a particular band of interest served as 
a starting point. This set was reduced immediately 
to those transitions whose frequencies were within 
10 GHz (the approximate combined bandpass limit 
of the HgCdTe detector-rf-amplifier that was used) 
of the frequency of a CO2 laser or a CO laser tran- 
sition frequency. Those transitions which were 
blended with nearby transitions from other bands 
or isotopic species of the molecule of interest were 
deleted from this subset. An attempt was then 
made to select (from the candidates available at 
this point) and measure those transitions which 
permitted the best determination of the constants. 
As a minimal set, these included low-/ transitions 
in both the P- and R-branches in order to deter- 
mine the band center, intermediate-/ transitions 
for determining the B -value, and high-/ (60 to 100 
depending on the particular band strength) transi- 
tions for the centrifugal distortion constants. If 
possible, additional measurements were made in 
order to cover the entire band with the smallest 
gaps possible. This served to increase the redun- 
dancy and to minimize extrapolated values in the 
calculated frequencies. 

In some instances, it was not possible to realize 
the above goals. The region of interest often cov- 
ered 100 cm -1 or so and typical TDL frequency 
coverage (guaranteed by the vendor) was 15 cm -1 , 
although many times the coverage turned out to be 
larger. Within the 15 cm -1 region, there were usu- 
ally holes in coverage and it was necessary to buy 
the TDLs in pairs. An additional factor which was 
not within the experimenter's control concerned 
frequency holes (regions where the beatnote was 
not discernible) in the bandpass of the combination 
of the HgCdTe detector and the rf amplifiers 
which follow it. These two factors were the final 
restrictions on the set of measurements which were 
used to determine the spectroscopic constants for 
the band. 
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2.2 CO2 Laser Standard Frequencies 

The 1.25 m long C0 2 lasers used in these mea- 
surements were constructed by the late F. R. 
Petersen [2.1]. These had gratings for line selection 
and output mirrors typically coated for 85% trans- 
mission. In order to render the output beams 
nearly parallel, compensated output coupling mir- 
rors were utilized. By compensated, we mean that 
the anti-reflection coated surface of the mirror was 
formed with a smaller radius of curvature than the 
reflecting concave surface. Irises at both ends were 
available for mode discrimination and control of 
the power output, which was 1 to 3 W. These lasers 
were equipped with internal absorption cells (filled 
to a pressure of 5.3 Pa (40 mTorr) of carbon diox- 
ide) to provide frequency stabilization by the 
Freed-Javan technique [2.2]. Although some seven 
isotopic combinations of CO2 have been the object 
of extensive frequency measurements, only three 
isotopes have been used for the NIST heterodyne 
frequency measurements. These were 12 C 16 2 
(626), 13 C 16 2 (636) and I2 C 18 2 (828). (The short- 
hand notation used here involves the last digits in 
the rounded nuclear masses for the three atoms, O, 
C, and O. For example, the 16 12 C 16 laser is desig- 
nated (626).) In addition to the 1.25 m lasers, a 2 m 
laser with a similar grating and compensated out- 
put mirror (but with an external absorption cell for 
stabilization) was used for high-/ and hot band 
transitions. This longer laser has been operated 
with the three different isotopic gases. 

The frequencies of the CO2 transitions were ini- 
tially related to the cesium standard in an NBS ex- 
periment published in 1973 (Evenson et al. [2.1]). 
This was an experiment in which the frequency of 
the methane stabilized HeNe laser was measured. 
A new wavelength measurement (Barger and Hall 
[2.3]) of the same methane transition was concur- 
rently completed. The combination of these two 
quantities led to a new (at that time) value for the 
speed of light, 299 792 456.2(1.1) m/s [2.4]. 

Figure la shows a block diagram of the fre- 
quency chain used in the methane frequency mea- 
surement. For the sake of brevity, the values of the 
offset frequencies in that experiment are not dis- 
cussed here; the interested reader is referred to the 
original paper [2.1]. The essential details of the 
chain are as follows. An X-band klystron was 
phase-locked to a quartz crystal oscillator and the 
X-band frequency was accurately measured by a 
counter which was referred to the cesium standard. 
The seventh harmonic of the X-band frequency 
was used to phase-lock a 74 GHz klystron. The 
twelfth harmonic of the 74 GHz oscillator was used 



as a reference to frequency lock an HCN laser. The 
twelfth harmonic of the HCN laser frequency was 
close enough to the frequency of the 28 fim water 
vapor laser for the difference to be made up by the 
frequency from a phase-locked klystron. Three har- 
monics of the water vapor laser and another phase- 
locked klystron frequency moved the chain up to 
the CO2 laser region, to the transition Ru(10). An- 
other frequency from a phase-locked klystron 
moved the chain from the 9.3 ju.m band to the 10.6 
ju.m band, to Ri(30) in particular. Three harmonics 
of the frequency of the Ri(30) transition brought 
the chain up near the frequency of the methane- 
stabilized HeNe laser. (The measured values for 
these transitions are indicated in the boxes in Fig. 
la and the values are given in THz.) Petersen et al. 
[2.5] subsequently used the measured values for 
these two transitions to make additional CO2 laser 
frequency measurements and to generate tables of 
transitions for both the I and II branches. 

Since 1973 several different international labora- 
tories have repeated this frequency measurement 
experiment with minor variations and confirmative 
measurements have led to a new value for the 
speed of light which is now defined to be 299 792 
458 m/s [2.6]. The proposal and expected accep- 
tance of this new definition provided impetus to 
extend frequency measurements to the visible por- 
tion of the spectrum. Major efforts at NIST- 
Boulder accomplished this objective and experi- 
ments were published in 1983 by Pollock et al. 
[2.7], and by Jennings et al. [2.8]. 

Figure lb shows a diagram of the chain used to 
measure the iodine transition at 520 THz (576 nm) 
[2.7]. The details of this chain are of secondary im- 
portance for purposes here, but a sketch will be 
provided. The starting point was the stabilized 
HeNe laser at 3.39 ju.m. A best value for this transi- 
tion was determined from the most accurate results 
of the international measurements. The Pi(50) 
transition of a stabilized I3 CC>2 laser was measured 
relative to the stabilized HeNe laser, which was as- 
sumed to oscillate at the frequency 88 376 
181.609 ±0.009 MHz. The frequency resettability 
of the CO2 laser was reported to be within 5 parts 
in 10", which is better than the 1 part in 10 10 uncer- 
tainty in the frequency of the HeNe laser. (Pe- 
tersen et al. [2.9] later used this accurately mea- 
sured Pi(50) transition to measure other C0 2 tran- 
sitions and to generate improved calibration ta- 
bles.) 

To continue with the sketch, we note that a 
second I3 CC>2 transition Pi (52) was phase-locked to 
Pi(50) by means of a stabilized 62 GHz klystron. 
The synthesis scheme used to arrive at the CCL 
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Diagrams of schemes relating C0 2 frequencies to the cesium stan- 



frequency is indicated in Fig. lb. The remainder of 
the chain was locked from the top down. That is, 
the second harmonic of the 1.15 jxm HeNe laser 
frequency was locked to the frequency of the dye 
laser, which was in turn locked to the frequency of 
the o hyperfine component of the visible (576 nm) 
127 I 2 17-1 P(62) transition. The CCL was in turn 
locked to the frequency of the 1.15 fim HeNe 
laser. Both the HeNe and the CCL served as trans- 
fer oscillators in this scheme. The comparison was 
made at the CCL-CO2 laser point of the chain. In a 
separate experiment, the HeNe laser was locked to 
its Lamb dip and that frequency was determined. 
The 13 C0 2 frequency listed in Fig. lb is the fre- 
quency used in the 520 THz determination (it was 
frequency-offset-locked from a stabilized CO2 laser 
to remove the dither) and is not the frequency of 
the center of the transition. The values for the two 
higher frequencies are given in Fig. lb. Additional 
details may be found in Ref. [2.7], 



The frequencies currently used for the 12 C I6 02 
isotope are based on the most recent values (which 
made use of the above results) which were pub- 
lished in 1983 by Petersen et al. [2.9]. The stated 
1 a uncertainties for the calculated tables based on 
these measurements are smaller than 5 kHz for the 
CO2 transitions (/ < 40) which were used in the in- 
frared heterodyne measurements. Subsequent CO2 
measurements relative to the values of Petersen 
were made by Freed and coworkers at MIT [2.10]. 
The uncertainties in the MIT values are also less 
than 5 kHz for the 13 C 16 C>2 transitions (and less 
than 10 kHz for the l2 C l8 02 transitions) used here. 
Note that the 1.25 m C0 2 lasers used in the TDL 
measurements were first generation lasers and that 
the numbers given in Tables 1, 2, and 3 were ob- 
tained using second generation lasers. Even if the 
reproducibility in the laser lock for the older lasers 
were somewhat poorer (an uncertainty of about 50 
kHz was allowed for the realization of the frequen- 
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cies) it would be of negligible consequence for the 
results presented here. The transition frequency 
values (for the three carbon dioxide isotopes) 
which were used are given in Tables 1, 2, and 3. 

Listed in Table 4 are the frequencies and 
wavenumbers for some transitions of the 13 C02 hot 



band (Ol'HllWOl) which extend the useful 
range of standard reference frequencies. These fre- 
quencies were determined from NIST measure- 
ments in which the 2 m CO2 laser and an external 
reference cell were used. While measurements 
were made on both the 636 and 626 isotopes [2.11, 



Table 1. Frequencies (MHz) for the 626 carbon dioxide laser 



Table 2. Frequencies (MHz) for the 636 carbon dioxide laser 



Rot. 


Frequency 


Rot. 


Frequency 


Rot. 


Frequency 


Rot. 


Frequency 


trans. 


(MHz) 


trans. 


(MHz) 


trans. 


(MHz) 


trans. 


(MHz) 


P(50) 


27 413 600.4112 


P(50) 


30 480 527 .0432 


P(50) 


26 035 339.9907 


P(50) 


29 076 007.9206 


P(48) 


27 478 430.1483 


P(48) 


30 545 874.3287 


P(48) 


26 096 450 .6641 


P(48) 


29 143 127.3150 


P(46) 


27 542 482.6299 


P(46) 


30 610 516.1385 


P(46) 


26 156 946 .4184 


P(46) 


29 209 472.6217 


P(44) 


27 605 762.5803 


P(44) 


30 674 445.7649 


P(44) 


26 216 830 .6112 


P(44) 


29 275 036.8793 


P(42) 


27 668 274.4491 


P(42) 


30 737 656.7009 


P(42) 


26 276 106 .3711 


P(42) 


29 339 813.3299 


P(40) 


27 730 022.4165 


P(40) 


30 800 142.6462 


P(40) 


26 334 776 .6055 


P(40) 


29 403 795.4269 


P(38) 


27 791 010.3989 


P(38) 


30 861 897.5131 


P(38) 


26 392 844 .0080 


P(38) 


29 466 976.8408 


P(36) 


27 851 242.0547 


P(36) 


30 922 915.4319 


P(36) 


26 450 311 .0648 


P(36) 


29 529 351.4663 


P(34) 


27 910 720.7882 


P(34) 


30 983 190.7566 


P(34) 


26 507 180 .0614 


P(34) 


29 590 913.4283 


P(32) 


27 969 449.7554 


P(32) 


31 042 718.0701 


P(32) 


26 563 453 .0887 


P(32) 


29 651 657.0881 


P(30) 


28 027 431.8676 


P(30) 


31 101 492.1893 


P(30) 


26 619 132 .0481 


P(30) 


29 711 577.0488 


P(28) 


28 084 669.7958 


P(28) 


31 159 508.1695 


P(28) 


26 674 218 .6570 


P(28) 


29 770 668.1612 


P(26) 


28 141 165.9746 


P(26) 


31 216 761.3094 


P(26) 


26 728 714 .4536 


P(26) 


29 828 925.5288 


P(24) 


28 196 922.6060 


P(24) 


31 273 247.1550 


P(24) 


26 782 620.8011 


P(24) 


29 886 344.5126 


P(22) 


28 251 941.6621 


P(22) 


31 328 961.5037 


P(22) 


26 835 938 .8920 


P(22) 


29 942 920.7360 


P(20) 


28 306 224.8891 


P(20) 


31 383 900.4083 


P(20) 


26 888 669 .7515 


P(20) 


29 998 650.0890 


P(18) 


28 359 773.8096 


P(18) 


31 438 060.1801 


P(18) 


26 940 814 .2413 


P(18) 


30 053 528.7321 


P(16) 


28 412 589.7252 


P(16) 


31 491 437.3923 


P(16) 


26 992 373 .0624 


P(16) 


30 107 553.1003 


P(14) 


28 464 673.7190 


P(14) 


31 544 028.8828 


P(14) 


27 043 346 .7579 


P(14) 


30 160 719.9062 


P(12) 


28 516 026.6578 


P(12) 


31 595 831.7569 


P(12) 


27 093 735 .7156 


P(12) 


30 213 026.1432 


P(10) 


28 566 649.1936 


P(10) 


31 646 843.3897 


P(10) 


27 143 540.1699 


P(10) 


30 264 469.0880 


P(8) 


28 616 541.7658 


P(8) 


31 697 061.4282 


P(8) 


27 192 760.2040 


P(8) 


30 315 046.3033 


P(6) 


28 665 704.6019 


P(6) 


31 746 483.7926 


P(6) 


27 241 395.7512 


P(6) 


30 364 755.6398 


P(4) 


28 714 137.7193 


P(4) 


31 795 108.6785 


P(4) 


27 289 446.5964 


P(4) 


30 413 595.2373 


P(2) 


28 761 840.9258 


P(2) 


31 842 934.5572 


P(2) 


27 336 912.3769 


P(2) 


30 461 563.5271 


R(0) 


28 832 026.2179 


R(0) 


31 913 172.5750 


R(0) 


27 407 012.8973 


R(0) 


30 531 879.5457 


R(2) 


28 877 902.4362 


R(2) 


31 958 996.0676 


R(2) 


27 453 013.4681 


R(2) 


30 577 664.6182 


R(4) 


28 923 046.4283 


R(4) 


32 004 0173874 


R(4) 


27 498 426.5523 


R(4) 


30 622 575.1932 


R(6) 


28 967 457.0638 


R(6) 


32 048 236.2545 


R(6) 


27 543 251.1292 


R(6) 


30 666 611.0177 


R(8) 


29 011 133.0037 


R(8) 


32 091 652.6661 


R(8) 


27 587 486.0315 


R(8) 


30 709 772.1308 


R(10) 


29 054 07Z6995 


R(10) 


32 134 266.8957 


R(10) 


27 631 129.9443 


R(10) 


30 752 058.8623 


R(12) 


29 096 274.3924 


R(12) 


32 176 079.4916 


R(12) 


27 674 181.4045 


R(12) 


30 793 471.8321 


R(14) 


29 137 736.1122 


R(14) 


32 217 091.2759 


R(14) 


27 716 638.7993 


R(14) 


30 834 011.9476 


R(16) 


29 178 455.6756 


R(16) 


32 257 303.3427 


R(16) 


27 758 500.3646 


R(16) 


30 873 680.4025 


R(18) 


29 218 430.6853 


R(18) 


32 296 717.0558 


R(18) 


27 799 764.1833 


R(18) 


30 912 478.6740 


R(20) 


29 257 658.5273 


R(20) 


32 335 334.0465 


R(20) 


27 840 428.1829 


R(20) 


30 950 408.5203 


R(22) 


29 296 136.3697 


R(22) 


32 373 156.2114 


R(22) 


27 880 490.1333 


R(22) 


30 987 471.9773 


R(24) 


29 333 861.1596 


R(24) 


32 410 185.7086 


R(24) 


27 919 947.6441 


R(24) 


31 023 671.3556 


R(26) 


29 370 829.6209 


R(26) 


32 446 424.9556 


R(26) 


27 958 798.1612 


R(26) 


31 059 009.2364 


R(28) 


29 407 038.2514 


R(28) 


32 481 876.6251 


R(28) 


27 997 038.9638 


R(28) 


31 093 488.4680 


R(30) 


29 442 483.3197 


R(30) 


32 516 543.6414 


R(30) 


28 034 667.1602 


R(30) 


31 127 112.1610 


R(32) 


29 477 160.8619 


R(32) 


32 550 429.1766 


R(32) 


28 071 679.6844 


R(32) 


31 159 883.6838 


R(34) 


29 511 066.6779 


R(34) 


32 583 536.6463 


R(34) 


28 108 073.2910 


R(34) 


31 191 806.6579 


R(36) 


29 544 196.3277 


R(36) 


32 615 869.7049 


R(36) 


28 143 844.5509 


R(36) 


31 222 884.9524 


R(38) 


29 576 545.1272 


R(38) 


32 647 432.2414 


R(38) 


28 178 989.8459 


R(38) 


31 253 122.6787 


R(40) 


29 608 108.1437 


R(40) 


32 678 228.3735 


R(40) 


28 213 505.3631 


R(40) 


31 282 524.1845 


R(42) 


29 638 880.1914 


R(42) 


32 708 262.4432 


R(42) 


28 247 387.0891 


R(42) 


31 311 094.0480 


R(44) 


29 668 855.8266 


R(44) 


32 737 539.0112 


R(44) 


28 280 630.8034 


R(44) 


31 338 837.0715 


R(46) 


29 698 029.3421 


R(46) 


32 766 062.8508 


R(46) 


28 313 232.0718 


R(46) 


31 365 758.2751 


R(48) 


29 726 394.7621 


R(48) 


32 793 838.9426 


R(48) 


28 345 186.2387 


R(48) 


31 391 862.8896 


R(50) 


29 753 945.8362 


R(50) 


32 820 872.4682 


R(50) 


28 376 488.4197 


R(50) 


31 417 156.3497 
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2.12], only the 636 frequencies were used for subse- 
quent TDL measurements. The 2 cr uncertainties 
for the transitions used for OCS measurements 
were less than 80 kHz. 

For the measurements in the CO2 laser region, 
the CO2 frequencies given in some of our papers 

Table 3. Frequencies (MHz) for the 828 carbon dioxide laser 



were rounded to the nearest 0.1 MHz. However, 
the full accuracy of the CO2 frequencies was re- 
tained for calculating or synthesizing the CO fre- 
quencies and the CO frequency was then rounded 
to the nearest 0.1 MHz. 



Rot. 


Frequency 


Rot. 


Frequency 


trans. 


(MHz) 


trans. 


(MHz) 


P(50) 


27 702 788.3340 


P(50) 


31 275 461.0484 


P(48) 


27 763 916.5271 


P(48) 


31 330 648.1933 


P(46) 


27 824 219.8818 


P(46) 


31 385 304.7512 


P(44) 


27 883 701.5085 


P(44) 


31 439 427.2029 


P(42) 


27 942 364.3439 


P(42) 


31 493 012.1188 


P(40) 


28 000 211.1532 


P(40) 


31 546 056.1633 


P(38) 


28 057 244.5316 


P(38) 


31 598 556.0988 


P(36) 


28 113 466.9070 


P(36) 


31 650 508.7885 


P(34) 


28 168 880.5416 


P(34) 


31 701 911.2012 


P(32) 


28 223 487.5337 


P(32) 


31 752 760.4139 


P(30) 


28 277 289.8196 


P(30) 


31 803 053.6153 


P(28) 


28 330 289.1753 


P(28) 


31 852 788.1092 


P(26) 


28 382 487.2180 


P(26) 


31 901 961.3173 


P(24) 


28 433 885.4075 


P(24) 


31 950 570.7820 


P(22) 


28 484 485.0477 


P(22) 


31 998 614.1692 


P(20) 


28 534 287.2879 


P(20) 


32 046 089.2707 


P(18) 


28 583 293.1240 


P(18) 


32 092 994.0069 


P(16) 


28 631 503.3995 


P(16) 


32 139 326.4282 


P(14) 


28 678 918.8066 


P(14) 


32 185 084.7176 


P(12) 


28 725 539.8870 


P(12) 


32 230 267.1924 


P(10) 


28 771 367.0327 


P(10) 


32 274 872.3053 


P(8) 


28 816 400.4870 


P(8) 


32 318 898.6464 


P(6) 


28 860 640.3445 


P(6) 


32 362 344.9440 


P(4) 


28 904 086.5522 


P(4) 


32 405 210.0657 


P(2) 


28 946 738.9095 


P(2) 


32 447 493.0190 


R(0) 


29 009 228.1753 


R(0) 


32 509 824.0588 


R(2) 


29 049 894.0639 


R(2) 


32 550 648.1734 


R(4) 


29 089 764.2422 


R(4) 


32 590 887.7557 


R(6) 


29 128 837.8481 


R(6) 


32 630 542.4476 


R(8) 


29 167 113.8723 


R(8) 


32 669 612.0318 


R(10) 


29 204 591.1583 


R(10) 


32 708 096.4309 


R(12) 


29 241 268.4016 


R(12) 


32 745 995.7070 


R(14) 


29 277 144.1494 


R(14) 


32 783 310.0604 


R(16) 


29 312 216.8002 


R(16) 


32 820 039.8289 


R(18) 


29 346 484.6028 


R(18) 


32 856 185.4857 


R(20) 


29 379 945.6557 


R(20) 


32 891 747.6385 


R(22) 


29 412 597.9061 


R(22) 


32 926 727.0276 


R(24) 


29 444 439.1492 


R(24) 


32 961 124.5238 


R(26) 


29 475 467.0268 


R(26) 


32 994 941.1261 


R(28) 


29 505 679.0261 


R(28) 


33 028 177.9601 


R(30) 


29 535 072.4788 


R(30) 


33 060 836.2746 


R(32) 


29 563 644.5594 


R(32) 


33 092 917.4396 


R(34) 


29 591 392.2837 


R(34) 


33 124 422.9433 


R(36) 


29 618 312.5075 


R(36) 


33 155 354.3890 


R(38) 


29 644 401.9248 


R(38) 


33 185 713.4920 


R(40) 


29 669 657.0662 


R(40) 


33 215 502.0763 


R(42) 


29 694 074.2966 


R(42) 


33 244 722.0715 


R(44) 


29 717 649.8141 


R(44) 


33 273 375.5085 


R(46) 


29 740 379.6471 


R(46) 


33 301 464.5166 


R(48) 


29 762 259.6531 


R(48) 


33 328 991.3192 


R(50) 


29 783 285.5157 


R(50) 


33 355 958.2301 



Table 4. Frequencies and 
-[11 'O.OS'O], band of 13 C0 2 



wavenumbers for the 01 '1 



Rot. 
trans. 



Frequency 
(MHz) 



P(50) 25 

P(49) 25 

P(48) 25 

P(47) 25 

P(46) 25 

P(45) 25 

P(44) 25 

P(43) 25 

P(42) 25 

P(41) 25 

P(40) 25 

P(39) 25 

P(38) 25 

P(37) 25 

P(36) 25 

P(35) 25 

P(34) 25 

P(33) 25 

P(32) 25 

P(31) 25 

P(30) 25 

P(29) 25 

P(28) 25 

P(27) 25 

P(26) 25 

P(25) 25 

P(24) 25 

P(23) 25 

P(22) 25 

P(21) 25 

P(20) 25 

P(19) 26 

P(18) 26 

P(17) 26 

P(16) 26 

P(15) 26 

P(14) 26 

P(13) 26 

P(12) 26 

P(ll) 26 

P(10) 26 

P(9) 26 

P(8) 26 

P(7) 26 

P(6) 26 

P(5) 26 

P(4) 26 

P(3) 26 

P(2) 26 



110 914.276(1909) 
163 452.505(1658) 
173 184.219(1500) 
223 707.197(1300) 
234 801.135(1157) 

283 388.000(1001) 
295 767.670(874) 
342 496.069(755) 

356 086.324(643) 
401 032.476(555) 
415 759.456(459) 
458 998.211(395) 
474 789.277(315) 
516 394.177(271) 
533 177.855(206) 
573 221.199(177) 
590 927.113(127) 
629 480.014(109) 
648 038.830(72) 
685 171.278(62) 
704 514.640(39) 
740 295.564(33) 
760 356.034(24) 
794 853.362(22) 
815 564.357(21) 
848 845.075(21) 
870 140.808(21) 
902 271.029(20) 
924 086.446(20) 
955 131.461(18) 
977 402.180(20) 
007 426.528(18) 
030 088.780(21) 
059 156.302(20) 
082 146.866(21) 
110 320.773(28) 
133 576.918(20) 
160 919.848(39) 
184 379.270(22) 
210 953.348(55) 
234 554.110(33) 
260 421.015(76) 

284 101.484(57) 
309 322.503(103) 
333 021.291(91) 

357 657.387(136) 
381 313.287(134) 
405 425.156(176) 
428 977.084(187) 



Rot. 

trans. 



Frequency 
(MHz) 8 



R(l) 26 

R(2) 26 

R(3) 26 

R(4) 26 

R(5) 26 

R(6) 26 

R(7) 26 

R(8) 26 

R(9) 26 

R(10 

R(ll 

R(12 

R(13 

R(14 

R(15 

R(16 

R(17 

R(18 

R(19 

R(20 

R(21 

R(22 

R(23' 

R(24 

R(25 

R(26; 

R(27 

R(28 

R(29 

R(30 

R(31 

R(32 

R(33 

R(34 

R(35 

R(36 

R(37 

R(38 

R(39; 

R(40; 

R(41 

R(42 

R(43' 

R(44 

R(45 

R(46 

R(47 

R(48 

R(49 

R(50; 



522 359.352(307) 
545 384.349(356) 
568 137.000(372) 
590 844.300(435) 
613 344.206(446) 
635 672.514(520) 
657 980.003(528) 
679 867.807(608) 
702 043.339(622) 
723 428.851(700) 
745 533.080(728) 
766 354.173(792) 
788 448.007(850) 
808 642.157(885) 
830 786.820(992) 
850 291.041(979) 
872 548.133(1159) 
891 298.919(1072) 
913 730.480(1356) 
931 663.740(1167) 
954 332.308(1589) 
971 383.309(1267) 
994 351.983(1865) 
010 455.287(1377) 
033 787.788(2192) 
048 877.189(1503) 
072 637.922(2577) 
086 646.387(1656) 
110 900.501(3027) 
123 760.107(1847) 
148 573.557(3552) 
160 215.432(2089) 
185 655.039(4159) 
196 009.301(2394) 
222 142.814(4856) 
231 138.505(2775) 
258 034.663(5653) 
265 599.695(3242) 
293 328.288(6559) 
299 389.374(3806) 
328 021.302(7584) 
332 503.902(4475) 
362 111.241(8736) 
364 939.495(5260) 

395 595.552(****) 

396 692.223(6167) 
428 471.603(****) 
427 758.013(7208) 
460 736.676(***») 
458 132.647(8390) 



" The number in parentheses is the estimated 1 a uncertainty in 
the last digits. 
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23 Heterodyne Frequency Measurements with 
the TDL and CO* Laser (860 to 1120 cm 1 ) 

The measurement procedure has evolved over 
the course of this work [5.73, 5.83, 5.87, 5.221, 
5.230, 5.243]. A brief history of the TDL refrigera- 
tion evolution and the details of the apparatus and 
procedure currently in use will be described here. 
For considerations involving the TDL, refer to 
Fig. 2, which is a block diagram of the measure- 



ment scheme recently used for some N2O measure- 
ments [5.243]. 

The first commercially available TDL spectrom- 
eters featured liquid helium Dewars as the refrig- 
eration system. Our initial system used a 4 L 
helium Dewar. It was necessary to use an assort- 
ment of stainless steel shims between the Dewar's 
OFHC (oxygen free high conductivity) copper cold 
surface and the TDL in order to vary the tempera- 
ture (and operating wavenumber) of the TDL. The 
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use of liquid helium was inconvenient, changing 
shims was cumbersome, and the resulting tempera- 
ture cycling of the TDLs was reputed to shorten 
their lifetime; nevertheless a narrow TDL 
linewidth was generally observed. 

As the demand for TDLs increased, vendors be- 
gan experiencing difficulties in growing semicon- 
ductors that would meet the customers' specified 
frequency region while operating in the 4 to 10 K 
range. The materials problem became more 
tractable as the temperature constraint was re- 
moved by selling to the customers closed-cycle 
coolers with 10 to 70 K operating capabilities. Soon 
it became nearly impossible to obtain TDLs opera- 
ble at helium temperatures. Most of our measure- 
ments were made using a closed-cycle cooler. This 
was an improvement in many areas, however in 
spite of an isolation scheme, some residual vibra- 
tions from the cooler's piston were transmitted to 
the TDL, resulting in the famous jitter linewidth 
which is familiar to all TDL users. The vendor has 
made two additional isolation improvement 
schemes available and while these appear worth- 
while, they are still not the ultimate solution. 

Currently, the NIST liquid helium Dewar has 
been modified to accommodate a four-laser mount- 
ing platform (including the heater coils and tem- 
perature sensing diodes) from a closed-cycle 
cooler. While helium consumption is too high for 
most of the TDLs available to us, this does offer 
promise with the recent advent of the higher tem- 
perature MBE (molecular beam epitaxy) TDLs. A 
few of these are now available and operate in the 
temperature range accessible with liquid nitrogen 
(70 to 120 K if a heater is available). 

Both the Dewar and closed-cycle refrigerator are 
interchangeable in that they were compatible with 
the laser control module (current controller) and 
temperature control system. Both the control mod- 
ule and temperature controller have also been up- 
graded to reduce current noise and temperature 
instabilities; these upgrades have proved worth- 
while. The particular type of refrigeration is not 
specified in Fig. 2. 

After passing through the AR-coated ZnSe win- 
dow of the refrigeration stage, the TDL radiation 
was collimated with an AR-coated //I lens and di- 
rected with a flat mirror and off-axis parabolic mir- 
ror into a 0.8 m Ebert-Fastie monochromator. A 
second off-axis parabolic mirror recollimated the 
TDL beam after it emerged from the monochroma- 
tor. A portion of the beam was split off and passed 
through an absorption cell (containing the 
molecule of interest, N 2 in this case) to a detector 



which was used initially in recording the spectra 
and later for the TDL locking procedure. The 
monochromator and a solid 3 in germanium etalon 
were used to help identify the particular molecular 
transition of interest. Once the transition had been 
identified, the kinematic mirror mount directing 
the TDL beam through the etalon was removed 
and the TDL beam was then focused on the 
HgCdTe mixer element. At this time, both the en- 
trance and exit slits were removed from the 
monochromator, in order to eliminate the fringes 
or channel spectra which the slits may cause via 
feedback to the laser. 

At this point in the procedure, it was useful to 
measure the TDL linewidth by heterodyning its 
output with that from a CO2 laser (or CO laser as 
described in the next section). TDL linewidths of 
several hundred megahertz are not uncommon at 
higher currents (higher gains) when viewed for in- 
tegration times of several seconds. Next, the cur- 
rent was reduced (while increasing the tempera- 
ture to maintain frequency) until an acceptable (10 
to 30 MHz), or at least the narrowest attainable, 
linewidth was achieved for a measurement. 



2.4 TDL Locking Procedure and Minimization of 
Error 

Assuming the best conditions, that is, a strong 
TDL (0.5 mW for example) with a single mode that 
has a flat power versus wavelength curve, a narrow 
TDL linewidth, and a well isolated and intense line 
for a locking reference, we could likely make a 
measurement with an uncertainty of less than 1 
MHz. While these conditions sometimes prevail, 
more often, they do not. 

It is germane to discuss here the TDL locking 
procedure and some of the ways we have mini- 
mized errors that can creep into a measurement. 
After the chopper was removed from the TDL 
beam path, the TDL was frequency-modulated at 

4.5 kHz and a first derivative scheme was used to 
lock the TDL frequency. In this procedure, the fre- 
quency of the radiation output of the TDL was 
tuned from well below the molecular transition to 
well above it, and the resulting derivative signal 
(including the baseline and the absorption line) 
was traced by the recorder. (See derivative trace in 
the xyy recorder box in Fig. 2.) The frequency was 
returned to the low value and then the recorder 
signal monitored as the trace was followed to the 
derivative midpoint for locking. (Note: The deriva- 
tive trace in Fig. 2 has a line labeled error signal 
superimposed on it. This error signal was recorded 
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while scanning the TDL current both up and down 
while the laser was locked with moderate gain. For 
infinite gain the line becomes horizontal, or zero 
volts everywhere; for nearly zero gain, the line has 
a slope near that of the derivative signal. The slope 
of the error signal may be chosen (by adjusting the 
loop gain) to any value between these limits. By 
proper adjustments of the gain value and the TDL 
current the TDL may be locked (or stabilized) to 
nonzero values.) A low gain was used in the lock 
loop and the recorder pen position was monitored 
during the measurement. One of the reasons for 
this was that the TDL mode generally was not flat 
over the region of the line of interest. Rather than 
locking the TDL to the zero-voltage point, the 
TDL frequency was locked to a point where the 
derivative signal crossed the existing baseline 
derivative. The sources of error to be avoided are 
not only sloping background (with resultant zero 
offset in the derivative signal), but also possible in- 
strumental zero offsets from the lock-in amplifier. 
Both sources become magnified when dealing with 
weaker TDL modes or low level absorption and a 
very high sensitivity on the lock-in amplifier, and 
the procedure outlined above was essential. 

Since the background slopes may have either sign 
and compensation for zero offsets may be either too 
large or too small, these errors are random when 
spread over measurements of several different 
lines. In the case of strong absorption lines and a 
powerful TDL mode, the difference between the 
lock point and V was generally negligible com- 
pared to other sources of error. 

Another obvious source of error in determining 
the location of the center of the molecular absorp- 
tion line was noise. In addition to TDL amplitude 
fluctuations and detector noise, feedback fringes or 
channel spectra are also included, although some of 
amplitude noise can be attributable to feedback. 
Several techniques were used to minimize sources 
of error. The frequency of the TDL modulation was 
chosen to be higher than that of the TDL amplitude 
fluctuations and higher than the upper frequency of 
typical detector noise. The usual techniques for 
fringe reduction, including the monochromator slit 
removal alluded to earlier, tilting of various optical 
elements in the TDL beam path (particularly the 
detectors) were all employed. Cells of absorbing 
gases (for isolation) were placed in the TDL beam 
path in a few instances although it was not possible 
for most spectral regions. 



2.5 Measurement of the Difference Frequency Be- 
tween the TDL and Gas (C0 2 or CO) Laser 

The portion of the TDL beam passing through 
the first beam splitter in Fig. 2 has its polarization 
rotated by 90° to be parallel to the gas laser beam 
polarization. (The TDL polarization was assumed 
to be in the plane of the junction; however, a size- 
able perpendicular component may also exist.) The 
TDL beam was then focused (typically, FL=12.5 
cm) through a second beam splitter onto the 
HgCdTe fast detector or mixer which had an ele- 
ment with an area of 0.1 mm 2 . The 3 dB bandwidth 
of this detector was 1 GHz. Power from the gas 
laser was focused with a 40 cm focal length lens and 
then reflected off the beam splitter (NaCl was cho- 
sen to keep the local oscillator power below 10 mW) 
in such manner as to make the gas laser beam 
col linear with the TDL beam and to make the beam 
waists coincide. 

After initial observation of the beatnote, the am- 
plitude was maximized by fine adjustment of the fo- 
cusing lenses. It was also necessary to ascertain that 
the beatnote observed on the spectrum analyzer 
was the one of interest. This was particularly rele- 
vant when using the CO laser. On some occasions, 
the first beatnote observed was due to the TDL ra- 
diation mixing with that from a nearby unintended 
CO transition which could not be prevented from 
lasing along with the CO transition of interest. Mul- 
timode TDLs are also the source of extraneous 
beatnotes. 

Once we determined that the observed beatnote 
was the one of interest, the TDL frequency was 
scanned by changing the current and the beatnote 
was followed on the spectrum analyzer from zero 
frequency up (or down) to the molecular feature to 
be measured. The scan rate was reduced to a lower 
value and the progress of the beatnote carefully 
monitored relative to the derivative signal. The 
TDL frequency was then locked to the desired 
point on the derivative signal which was displayed 
on the recorder. The beatnote was then averaged 
with the persistent screen averaging feature of the 
spectrum analyzer, and a marker oscillator was ad- 
justed to the center of this averaged display. (A rep- 
resentative beat note is shown in the right hand 
portion of Fig. 2. Here the frequency span of the 
spectrum analyzer display was 100 MHz.) The 
marker oscillator frequency was counted and the 
measurement repeated a number of times (10 to 20, 
depending on the reproducibility of the measure- 
ments). 
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2.6 Minimization of the Difference Frequency 
Uncertainties 

Ideally, the best method to determine the differ- 
ence frequency between the TDL and the gas laser 
would be to use an electronic counter. However, 
the signal-to-noise ratio (S/N) of the beatnote and 
the frequency modulation associated with the cold 
head and compressor generally precluded this ap- 
proach. The next best approach is that described in 
the preceding paragraph. 

The best measurements were those made with a 
liquid helium Dewar; the beatnote was essentially 
stationary on the spectrum analyzer. However, the 
rapid He consumption rates for higher tempera- 
ture TDLs made this choice impractical as well as 
inconvenient. When the compressor was used, the 
beatnote had a jitter linewidth associated with it 
and its frequency fluctuations made determination 
of the beatnote center more difficult. A wide varia- 
tion of jitter linewidths from many different TDLs 
has been observed over an extended period. The 
current tuning rates (60 to 1600 MHz/ma) and 
linewidths due to current noise vary widely from 
one TDL to another. In a similar fashion the jitter 
linewidth varies greatly from one laser to the next, 
due to varying sensitivity to vibrations associated 
with the compressor/coldhead. Sometimes an ap- 
parent jitter linewidth was due to feedback, how- 
ever this was generally recognized and steps were 
taken to minimize it. Beatnotes ranging in width 
from a few megahertz (10 ju,m TDLs in a liquid 
helium Dewar) to 60 to 100 MHz (6 ju.m TDLs in a 
closed-cycle cooler) were observed during the mea- 
surements. However the larger values (in the 5 to 6 
ju,m region) were observed prior to the currently 
implemented improvements in the vibration isola- 
tion system. The most recent approach was to ad- 
just the current modulation for the derivative lock 
such that the beatnote linewidth was not broad- 
ened beyond the jitter linewidth. (This was subject 
to retention of a suitable S/N for the lock signal.) 

The frequency-modulated beatnote was ob- 
served at a repetitive rate on the spectrum ana- 
lyzer. The pulse rate of the closed-cycle cooler was 
asynchronous with both the modulation rate and 
the spectrum analyzer sweep rate. As a result, the 
beatnote observed on the spectrum analyzer made 
small and slowly varying excursions about an aver- 
age value. This led to some scatter in the measured 
value for the beatnote center frequency. Some 
experiments were conducted to check the compres- 
sor-induced fluctuations as a source of systematic 
error. Typically, 20 measurements were made with 



the compressor on, and the marker oscillator fre- 
quency was adjusted to the resultant average value. 
The TDL lock point was then rechecked and the 
compressor was turned off (eliminating the jitter 
from the beatnote) momentarily. The jitter-free 
beatnote was observed for a few seconds in this 
configuration. To date no appreciable deviation of 
the jitter-free beatnote from the marker oscillator 
has been observed. 

On some occasions, the beatnote envelope was 
slightly asymmetric. Generally, two different opera- 
tors have determined the center value and some 
subjective disagreement was apparent. In a recent 
set of measurements of 20 transitions the average 
value of 10 measurements each from two operators 
varied by 2.5 MHz. This was well within the as- 
signed uncertainty of 7 MHz for the measurement. 
More often, the average values from different op- 
erators agree within a fraction of 1 MHz. 

Another difficulty in these measurements was 
the presence of holes in the frequency coverage of 
the detector-rf-amplifier combination. In some in- 
stances, these frequency holes were associated with 
connector lengths and their effect could be mini- 
mized (by moving the hole to another frequency) 
by using line stretchers or by changing cables. 
Some holes were associated with lengths of con- 
necting elements in the detector, and for practical 
purposes could not be eliminated. In other cases, 
holes were associated with amplifiers themselves; 
sometimes they precluded making measurements. 
In a few instances shallow holes have led to system- 
atic errors. This occurred when only very weak 
TDL modes were available and the S/N for the 
beatnote was small (3 to 4 dB for example). Often 
the beatnote envelope was fairly wide (50 MHz or 
greater). In cases like these, one side of the beat- 
note envelope can overlap a hole and the apparent 
line center will be shifted. This has happened in a 
few instances but the error was apparent in the fit- 
ting process. In these cases, a repetition of the 
measurement with a different TDL and a much 
stronger beatnote gave a different and better fitting 
result. Such holes generally remain at the same fre- 
quency and experience has shown which frequency 
regions to avoid. 

2.7 Measurements with a CO Laser Transfer 
Oscillator and C0 2 Laser Synthesizer 

Near the inception of this program, a CO laser 
stabilization scheme on low pressure CO laser dis- 
charges had been demonstrated by Freed [2.13]. 
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More recently, a stabilization scheme using opto- 
galvanic detection has been reported by Schneider 
et al. [2.14]; neither of these schemes was operable 
over the entire range required for our measure- 
ments. Some values of CO frequencies in the liter- 
ature (available when the measurements in this 
region began) were in error by over 50 MHz. Since 
the goal of this measurement program was to be 
able to make measurements with a 3 MHz uncer- 
tainty, it was necessary to measure the frequency of 
the CO laser at the same time that the CO laser- 
TDL difference frequency measurements were 
made on N2O and OCS transitions. This process 
required the use of the CO laser as a transfer oscil- 
lator, and the CO laser frequency was measured 
relative to a frequency generated by a CO2 laser 
synthesizer [5.94, 5.221, 2.15]. 

Two different CO lasers were used in this 
manner. One was a sealed-off laser which was 
cooled by flowing alcohol through dry ice and then 
through a jacket around the discharge tube. This 
operated over the frequency range from 1600 to 
1900 cm -1 (corresponding to lower vibrational 
quantum numbers ranging from w"=20 to v"=6). 
The second CO laser was a flowing gas laser which 
was cooled by liquid-nitrogen and operated from 
1220 to 1600 cm -1 (v =36 to v =20). After installa- 



tion of a shorter wavelength grating and an appro- 
priate output mirror, operation was extended to 
the 1900 to 2080 cm" 1 region (v =6 to v = 1). Ad- 
ditional details regarding the liquid-nitrogen 
cooled CO laser may be found in the literature 
[5.125, 5.224, 5.231]. 

Figure 3 shows a block diagram for making fre- 
quency measurements with the CO laser. The 
dashed outline shows the kinematically mounted 
mirror, Mi, in position for measuring the CO laser 
frequency relative to the CO2 synthesizer, which is 
shown enclosed in the large dashed box. The syn- 
thesizer consists of two stabilized CO2 lasers, a 
phase-locked microwave oscillator and frequency 
counter, a metal-insulator-metal (MIM) diode, and 
a combination of an rf amplifier, rf spectrum ana- 
lyzer, and a to 1.0 GHz rf frequency synthesizer. 
When radiation from the two CO2 lasers and the 
microwave oscillator were coupled to the MIM 
diode, currents were generated at a synthesized 
frequency, v s , given by 

where vi and m were the frequencies of the CO2 
laser frequency standards, and vu was a microwave 
frequency. The quantities /, m , and n are integers 
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which are allowed both positive and negative val- 
ues. The quantity (l + |/| + |m| + |«|) is called the 
mixing order; the synthesized currents generally 
become weaker as the mixing order is increased. 
Mixing orders vary from 3 or 4 near 50 THz to 7 or 
8 near 38 THz, the frequency at the longest wave- 
length operation of the CO laser used in these 
measurements. Typical values might be / = 3 or 4, 
m = — 2 or —3, and (with the use of an X-band 
klystron) n was restricted to 0, ±1, or ±2. 
When the CO laser radiation was focused on the 
MIM diode, an additional current at the CO laser 
frequency, i>to, was generated in the diode and it 
combined with the synthesized frequency, v s , to 
produce a difference frequency beatnote at a fre- 
quency, fm. (The microwave frequency was chosen 
so that this beatnote was within the 1.2 GHz range 
of the spectrum analyzer in use.) The beatnote was 
amplified, displayed on the spectrum analyzer, and 
its excursion was noted as the CO laser was tuned 
through its gain bandwidth. The beatnote was posi- 
tioned at the center of this excursion (a determina- 
tion of the frequency of the CO transition was a 
secondary objective) and a marker signal from the 
rf synthesizer was used to mark this frequency 
point on the spectrum analyzer. The rf synthesizer 
reading was then used as the value for i^i and the 
CO laser (which was not locked) was periodically 
readjusted to return the beatnote to the assigned 
frequency. 

The frequency of the CO laser (transfer oscilla- 
tor) was then 

^CO = transfer = V % ± Vb\. 

The full accuracy of CO2 frequencies was used for 
the v s calculation. An uncertainty of 0.3 MHz 
(which includes allowance for drift between read- 
justments) in the transfer oscillator was included in 
the measurement uncertainty of the molecular 
transition, which was given by 

Vmol= ^transfer± ^B2, 

where vbi was the beatnote between the TDL and 
the transfer oscillator. The transfer oscillator fre- 
quency and beatnote frequency are both rounded 
to the nearest 0.1 MHz. The main uncertainty was 
again due to the TDL linewidth which was dis- 
cussed in the early part of the chapter. 

The most recent advance in making heterodyne 
frequency measurements with TDLs involves a 
computer-controlled, frequency offset-locking (CC- 
FOL) scheme. Freed et al. [2.16] demonstrated the 
use of a frequency offset lock combined with a 



frequency synthesizer to control the output fre- 
quency of the TDL. We have combined that tech- 
nique with the scanning and data-logging 
technology used in this laboratory for other mea- 
surements [2.17-2.19] to obtain accurate data on 
dv/dP, the pressure-induced frequency shifts, in 
the rovibrational spectrum of OCS. The potential 
for better absolute frequency measurements was 
also demonstrated. 

Figure 4 shows a block diagram of the apparatus 
used for this type of measurement. The output 
beams from a TDL and a CO2 laser frequency stan- 
dard were focused with separate lenses and then 
combined with a ZnSe beam splitter and directed 
to a HgCdTe heterodyne mixer/detector which pro- 
duced a beatnote at the difference frequency, vb, 
between the two lasers. The beatnote was ampli- 
fied in an rf amplifier and displayed on an rf spec- 
trum analyzer. A balanced mixer was used to 
down-convert the beatnote at frequency vb to a 
nominal 160 MHz, the region of operation of the 
IF amplifier and discriminator. The beatnote was 
fed to one input arm of the balanced mixer, and 
the output of the sweepable frequency generator, 
at frequency v svl , was fed to the other input arm 
(the local oscillator arm). The frequency v svl was 
adjusted such that |hw-v b | was nominally 160 
MHz, and this resulting output signal was fed to 
the discriminator which had a sensitivity of 0.1 V7 
MHz and an 80 MHz bandwidth. 

After the switch in the loop filter was closed, the 
discriminator-based locking loop adjusted the TDL 
frequency to insure that the beatnote va was locked 
at a frequency »*>o a wa y from ^sw That is, 

pdo=\vsv/-vb\, 

where udo is close but (due to the presence of vari- 
ous zero offsets in the locking loop) not necessarily 
equal to 160 MHz. For frequency shift and 
lineshape measurements, the important point is 
that the frequency vdo must remain fixed, whatever 
value it assumes. If a frequency measurement is the 
objective, it becomes necessary to measure j^o. 
Frequency of the CCFOL TDL is then given by 

vtdl = ^co 2 ± Psw— vdo. 

Two frequency measurements were made with an 
uncertainty of ± 2 MHz, which was almost entirely 
due to the TDL linewidth. By narrowing the TDL 
linewidth with a faster loop filter, one should be 
able to use an electronic counter to measure v D o 
and make measurements with uncertainties the or- 
der of 0.2 MHz. 
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Fig. 4. Block diagram of a computer controlled frequency offset-locked spectrom- 
eter. 



We have restricted our initial experiments to 
those OCS transitions which lie within 2000 MHz 
of a CO2 laser transition because the lock loop 
requires a beatnote with a good SIN and the beat- 
note signal decreases with increasing frequency. 
This is also the band limit of our most convenient 
rf amplifier. We chose from the available TDLs 
those with sufficient power to give a beatnote with 
a SIN of about 30 dB. For the present measure- 
ments a 400 ms integration time was used and 640 
points were recorded in each direction. Recording 
in both directions is a good way to cancel certain 
types of systematic errors. Generally only one 
round trip pass was made per measurement but as 
many passes could be made as required to give a 
good SIN. 



A large number of measurements have been 
made with the older technique and these have been 
combined with FTS measurements, particularly the 
high quality measurements made recently. As the 
situation stands now good molecular constants ex- 
ist, with the band centers currently having the 
largest uncertainties. The number of measure- 
ments left to be made is a relatively small number 
of high quality. The best approach is to make sub- 
Doppler or saturated measurements, however, the 
accidental overlaps required are rather infrequent. 
We believe the CCFOL approach is the next best 
option and several strategic overlaps occur in the 2 
GHz range. 
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2.8 Measurements with the Color Center Laser 

Pollock et al. [5.219] used a color center laser 
(CCL) to perform a set of experiments on N 2 0. A 
brief description of their work and some related 
work concludes the summary of heterodyne tech- 
niques. It is of interest to compare and contrast 
some of the salient features of the TDL and the 
CCL. The tuning range of a TDL mode was 15 to 
30 GHz; that of the CCL was less than 1 GHz. The 
linewidth of the TDL in the best instances was a 
few MHz, that of the CCL was 10 kHz. Perhaps the 
most important feature of the CCL was a large 
power output (in excess of 10 mW), which along 
with its beam quality permitted a direct coupling of 
the CCL output to the MIM diode, and subsequent 
synthesis measurements without a transfer oscilla- 
tor. This relatively large power and narrow 
linewidth made it an ideal tool to use for some sat- 
urated absorption measurements by Pollock et al. 
on CO [5.304]. Sub-Doppler measurements could 
not be made on the N2O band studied, and uncer- 
tainties of 4 to 8 MHz (at 130 to 140 THz) were 
reported. This uncertainty was due to the uncer- 
tainty in locating the center of the transition and 
was due in part to the small free spectral range of 



the tuning element of the CCL, which frequently 
prevented us from sweeping over the entire line 
and was also insufficient to sweep far enough on 
either side of an absorption line to determine a 
background slope. 

Shown in Fig. 5 is a block diagram of the ring 
configuration color center laser developed by Pol- 
lock and Jennings [2.20]. The lasing entity in the 
lithium-doped potassium chloride crystal was an 
(F + 2)a center. The centers were optically pumped 
with the 3 W power output from an Nd:YAG laser 
operating in a TEMoo mode at 1.3 |xm. These color 
centers were continuously replenished by uv radia- 
tion from an Hg lamp. 

Two Brewster's angle sapphire prisms, a single 
plate birefringent filter and one etalon comprised 
the tuning elements. The ring was constrained to 
operate in a unidirectional manner by an optical 
diode consisting of an AR-coated YIG plate in a 
0.1 T magnetic field and a Brewster-cut quartz 
plate reciprocal rotator. One portion of the output 
radiation was used for stabilization to the side of a 
fringe in a passively stabilized optical cavity. The 
cavity was scanned by tuning a Galvo plate inside 
this reference cavity. Corrections were applied to 
the Galvo plate (slow) and to the PZT driving the 
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tuning mirror (fast) in the laser resonator; this nar- 
rowed the CCL linewidth to 10 kHz. A second por- 
tion of the beam was split off and sent through a 
cell containing CO. The Galvo plate in the refer- 
ence cavity was modulated at 7 kHz and slowly 
scanned to observe the first derivative signal, which 
was used to lock the CCL to the CO lines of inter- 
est. 

A third portion of the CCL radiation was di- 
rected to the CO2 synthesizer (more specifically the 
MIM diode portion) for a simultaneous measure- 
ment of the CO frequency. Typical synthesis 
schemes used 5^i, 4^i + vi, or 3v\ + vi, where v\ and 
vi are different CO2 laser frequencies. No mi- 
crowave oscillators were required, and the vb\ type 
beatnotes fell within 2 GHz. In contrast to the 
TDL measurements, the measurement uncertainty 
was entirely the uncertainty in locating the center 
of the absorption line. 
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3. Formulas and Data Sources Used to 
Prepare the Tables 

3.1 Expressions Used for Fitting the Frequency 
Data and for Calculating the Transition 
Wavenumbers 

For diatomic molecules and for linear tri- 
atomic molecules in 1 X electronic states the energy 
levels are generally given by 

E vJI =G v +B v J(J + l)-D v [J(J + l)-l 2 ] 2 

+//v[/(/ + l)-/ 2 ] J +L v [/(/ + l)-/ 2 ] 4 

+ higher terms, (3.1) 

where / is the quantum number for overall rota- 
tional angular momentum and / is the quantum 
number for vibrational angular momentum. Di- 
atomic molecules have no vibrational angular mo- 
mentum; that is, / = 0. In this work no higher order 
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terms were needed and even the H v and L v terms 
were either poorly determined or not deter- 
minable. 

For diatomic molecules, an alternative formula- 
tion is often given for the energy levels, 



£», = XY„(u + 1/2)' [J(/ + l)] / . 



(3.2) 



Dunham [3.1] has related the Y (; constants to the 
potential function of a diatomic molecule in a 1 X 
state. Most of the papers reporting constants for 
CO use Eq. (3.2). The ground state of the NO 
molecule is a 2 II state and is treated differently 
[5.371]. 

Transitional frequencies, Vcaic, are calculated as 
differences between energy levels 



Vcalc—E v jl'—E, 



yJl 



(3.3) 



In this book the band centers, v , are defined by 



vo — G v — G v ". 



(3-4) 



For linear triatomic molecules two types of per- 
turbations are commonly encountered that affect 
the importance of higher order terms in Eq. (3.1), 
/-type resonance and Fermi resonance. Because it 
has a large effect on the centrifugal distortion con- 
stants, /-type resonance was treated explicitly in 
the analysis of the OCS and N 2 spectra. Both /- 
type doubling and /-type resonance are manifesta- 
tions of the same matrix element that couples 
levels that differ only in the value of the / quantum 
number, where / is treated as a signed quantum 
number. In this book, both effects are treated un- 
der the general title of /-type resonance. If the 
bending vibrational quantum number, vi, is greater 
than zero, /-type resonance will be present and is 
usually noticeable. 

When V2 * 0, the /-type resonance was taken into 
account by diagonalizing the energy matrix which 
includes the matrix elements coupling / levels with 
/ ± 2 levels. The form of these matrices has been 
described in Refs. [5.120] and [3.2] but we shall 
repeat that description for a specific case. 

For U2 = 3, there are four possible values' of /, 
/ =3, / = 1, / = - 1, and / = -3. The /-doubling con- 
stant, q v , represented by 

?v=4v°-?vj/(/ + 1)+?vjj/ 2 (/ + 1) 2 , (3.5) 



couples these levels through the matrix element 



W,- 2 ,, = W,,i-2= <v,l\H\v,l -2> =l/4q v {[v 2 +l] 

[v 2 -I + 2]X[J(J + l)-l(l -1)][J(J + 1) 
-(/-l)(/-2)]}" 2 . (3.6) 



For each / level the form of the energy matrix for 
V2 = 3 is 



£°(v,/ = 3) 


W 3 ,i 


Wl.3 


E°(v,l = l) W t ,-i 





W-u E\v,l=-\) JF-1,-3 





W-i,-i £°(v,/=-3) 



(3.7) 



Here, the matrix elements are given by Eq. (3.1) 
(where E y ji =E°) and Eq. (3.6). / = is not allowed 
and for / < 3 only the central two-by-two matrix is 
allowed. Higher order terms coupling / and l±A 
levels are sometimes important but were not neces- 
sary for the present calculations. 

In general the /-type resonance calculation re- 
quires the use of a matrix of dimension vi + 1 by 
u 2 + l- Since this is a nonlinear system, a nonlinear 
least-squares fitting technique was needed to fit 
the experimental data to determine the best con- 
stants, as explained later. 

Most other workers have only used Eq. (3.1) to 
fit the data of OCS and N 2 0. That has the effect of 
absorbing the /-type resonance into the effective 
B v , D v , and H v values. While such a treatment is 
quite reasonable, the effective values of the higher 
order constants are quite different from the ground 
state values, and the level at which Eq. (3.1) is 
truncated has an important effect on both interpo- 
lation and extrapolation. By treating the /-type res- 
onance explicitly, we bring the effective values for 
Z?v and H v much closer to the ground state values. 
This model gives a better approximation to the true 
Hamiltonian than the model that only uses Eq. 
(3.1). This improvement in the model used for fit- 
ting the data improves the reliability of the least- 
squares fits and gives more accurate uncertainties 
for the calculated transition frequencies. 

For each value of |/| (except / =0) the states with 
V2 > are split into e and / components. For OCS 
and N 2 the / = states ( l X + states) always have 
the same symmetry (or parity) as the e states. 
These e and / components have been assigned in 
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accordance with the convention established by 
Brown et al. [3.3]. That convention leads to the se- 
lection rules: 

AJ=0, e+*f 

AJ=±1, e<^e, and/**/ 

for electric dipole transitions. These selection rules 
are obeyed even when the normal rule, Al = 0, ± 1, is 
broken because perturbations always connect e to e 
and / to /. 

All of the /-type resonance energy matrices, like 
Eq. (3.7), may be factored into two submatrices 
which represent the e levels in one case and the/ 
levels in the other case. We have used the full ma- 
trix, as indicated by Eq. (3.7), rather than a factored 
form, because it is more convenient for obtaining 
the eigenvectors needed to calculate the intensities 
of the transitions. 

The present analysis ignores the Fermi resonance 
that couples the levels (ui,i/2 + 2,/,f3 — 1,/) and 
(vi,V2,l, V3,J) of OCS and N 2 0. In OCS the unper- 
turbed Fermi resonance levels are far apart, so 
there is very little change in the resonance across a 
band. This results in only small changes in the effec- 
tive values of D y , H v , and L v . Such small changes can 
be accommodated by Eq. (3.1) without affecting ei- 
ther the accuracy of the least-squares fits or the ac- 
curacy of the calculated values. In N2O the Fermi 
resonance is expected to be more important but 
again the effective values of D v , H v , and L v are only 
slightly changed from the unperturbed values. 

Since the Fermi resonance coupling is different 
for different values of |/|, it gives rise to different 
effective values of the constants q v , B v , etc. for levels 
that differ only in the value of |/|. Consequently, in 
Eq. (3.7) one must use two values for By, one for the 
|/| = 1 states and one for the |/| = 3 states. Similarly, 
two values are needed for D y , H„, and L v . There will 
also be two different off-diagonal coupling con- 
stants, q v , in Eq. (3.7), one for the W^u W^, W-1,-3, 
and W-x-i terms and a slightly different one for the 
W\-\ and W-\,\ terms. However, these small differ- 
ences in q v are difficult to separate from the differ- 
ences in 5v and D». Consequently, we have been 
forced to use a single value of q v for a given vibra- 
tional state irrespective of differences in the value 
of the / quantum number. 

In analyzing the spectral data to get the rovibra- 
tional constants for calculating the most accurate 
transition frequencies of OCS and N2O, a large 
body of data on many different types of transitions 
was fit. Because of the form of the energy matrix for 



/-type resonance, it was not possible to use a linear 
least-squares technique to fit the data. Instead, an 
iterative nonlinear least-squares fitting procedure 
was used. In this procedure it was necessary to ap- 
proximate the derivative of the transition frequency 
with respect to each constant by applying the tech- 
nique described by Rowe and Wilson [3.4]. 

A similar nonlinear least-squares fitting proce- 
dure was used for NO but the energy matrix was 
somewhat different from Eq. (3.7). For a complete 
description of the energy matrix for NO one should 
refer to the work of Hinz et al. [5.371]. 

In order to calculate the statistical uncertainties 
in the calculated wavenumbers given in these tables 
it was necessary to use the variance-covariance ma- 
trix, given by the least-squares analysis, and the 
derivative of the transition frequency with respect 
to each constant. The uncertainty, or estimated 
standard error, given by a(v) was then determined 
by the double summation, 



<r(iO = S 2 Kij (dvldc)(dvldCj), 



(3.8) 



where Vy is a particular element of the variance- 
covariance matrix and dvldci and dv/dCj are the 
derivatives of the transition frequency with respect 
to the rovibrational constants c-, and q respectively. 

3.2 Data Sources Used for Fitting the Frequency 
Data and for Calculating the Transition 
Wavenumbers 

3.2.1 OCS All of the OCS transitions to be 
used for calibration (shown with an asterisk in the 
atlas) were calculated by means of constants and a 
variance-covariance matrix given by a single least- 
squares fit that included all of the frequency mea- 
surements given in the literature. The equations 
used in this fit were described in the preceding sec- 
tion. In this section we indicate what references 
provided the data that went into that fit and give a 
few more details about the fit and the selection of 
data. 

The rotational spectrum of OCS has been exten- 
sively studied by microwave and sub-millimeter 
wave techniques. These measurements use fre- 
quency techniques for calibration and have uncer- 
tainties on the order of ± 0.05 MHz and in some 
cases even smaller uncertainties. Such measure- 
ments are blessed with small line widths and are 
made at low pressures which contribute to the accu- 
racy of the measurements. The three most abun- 
dant isotopic species of OCS have no fine structure 
due to quadrupole effects. 
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Some microwave measurements [5.43, 5.79] ex- 
tend to fairly high / values so they are able to give 
accurate values for B y and D v . In addition, the het- 
erodyne measurements made by Vanek et al. 
[5.124] on high-/ transitions were used in the analy- 
sis. Although there are a great many measurements 
of rotational transitions for the lower vibrational 
states, Bogey and Bauer [5.78] and Tanaka et al. 
[5.98] have given measurements of rotational tran- 
sitions for fairly high vibrational states, up to 4100 
cm -1 . Some transitions show the splitting due to 
1-type resonance and for a few vibrational states, 
01*0 [5.32, 5.47], 02 2 [5.54], and 03*0 [5.54], the 
transitions between split levels have been observed. 

Altogether 333 frequency measurements of rota- 
tional transitions, taken from the above references 
as well as from Refs. [5.5, 5.29, 5.30, 5.37, 5.41, 
5.48, 5.55, 5.66, 5.119], were included in the least- 
squares fit that determined the rovibrational con- 
stants given in Tables 5 and 6. When possible these 
measurements were given uncertainties suggested 
in the original papers. In some cases the uncer- 
tainty was estimated by us, based on other work 
from that time or from that laboratory, or based on 
the goodness of the fit. 

With three exceptions, all of the infrared hetero- 
dyne frequency measurements came from a series 
of papers from the same laboratory at the National 
Institute of Standards and Technology in Boulder, 
Colorado [5.73, 5.83, 5.87, 5.94, 5.107, 5.120-5.122, 
5.125, 5.129]. The exceptions are the measure- 
ments on the 02°0-00°0 band and accompanying hot 
bands made in the Harry Diamond Laboratory by 
Sattler et al. [5.88], the measurements in the same 
frequency region made at the University of Lille by 
Fayt et al. [5.126], and the preliminary measure- 
ments from the University of Bonn [5.137a]. From 
the root-mean-square (rms) deviations of the NIST 
measurements it was obvious that their assigned 
uncertainties were too large by approximately a 
factor of two. This reflected caution in allowing for 
systematic errors which would not be revealed by 
the least-squares analysis. 

An extensive set of laser-Stark resonance mea- 
surements have been made by Fayt and others 
[5.105, 5.106, 5.111, 5.112] and by Tanaka et al. 
[5.71, 5.97]. This is a type of frequency measure- 
ment that should be quite accurate but we have not 
included those data in this least-squares fit as they 
would introduce the additional complications of 
determining dipole-moment functions and assess- 
ing the accuracy of the electric field measurements. 



In the case of measurements using CO laser transi- 
tions there is an additional uncertainty in the laser 
frequency. 

In order to determine the most accurate cen- 
trifugal distortion constants, some of the better 
diode laser and FTS measurements were included 
in the least-squares fits. For the most part the 
diode laser measurements given in Refs. [5.83, 
5.88, 5.94] were calibrated with the heterodyne 
measurements and probably had systematic errors 
much smaller than the dispersion shown by the 
least-squares fit. 

Except for the weakest transitions, the FTS mea- 
surements were more precise than any of the other 
infrared measurements. The high precision, how- 
ever, does not necessarily imply high accuracy. The 
FTS data used in the analysis for these tables were 
taken from Refs. [5.75, 5.101, 5.102, 5.120, 5.122, 
5.128, 5.132, 5.135, 5.136] or from private commu- 
nication with the authors of those papers in the 
cases where the original data were not published. 

The FTS measurements were given uncertainties 
equal to the rms deviations of the fits on a band-by- 
band basis. In order to keep the FTS measure- 
ments from affecting the determination of the 
vibrational energy levels, they were fitted to the 
same rotational constants as the other data, but to 
different band centers. All of the recommended 
calibration data were based on vibrational energy 
levels determined only from frequency measure- 
ments, not from FTS or ordinary diode-laser mea- 
surements. 

For the less abundant isotopomers of OCS many 
of the above papers plus a few additional papers 
[5.5, 5.52, 5.82, 5.84] give microwave and sub-mil- 
limeter wave measurements of rotational transi- 
tions. Some of the infrared heterodyne measure- 
ments also included transitions for the less abun- 
dant isotopomers of OCS [5.73, 5.87, 5.94, 5.120, 
5.122, 5.125]. A few transitions of 16 13 C 32 S and 
16 12 C 34 S have enough frequency measurements to 
warrant being considered as possible calibration 
transitions. In most cases, however, transitions of 
the rarer isotopomers are only included in the atlas 
to help in identifying the other, more useful transi- 
tions. 

A good many wavelength measurements have 
been made on the less abundant isotopomers, and 
they can be found in the bibliography, Sec. 5. The 
most important sources of information on infrared 
measurements of the less abundant species are 
Refs. [5.45, 5.65, 5.76, 5.101, 5.109, 5.111, 5.128]. 
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Table 5. 


Rovibrational constants (in cm -1 ) used for 16 12 C 32 S 








Vib. state 


Va 


fl. 


D.xlO" 


tf v xl0 M 


L.X10 17 


00"0 


0.0 


0.202 856 740 8(8)" 


4.340 64(25) 


-0.329(30) 




Ol'O 


520.422 055(147) 


0.203 209 834 8(21) 


4.411 48(31) 


-0.260(38) 




00"1 


858.966 932(48) 


0.202 251 831 6(60) 


4.433 50(36) 


0.045(45) 




02^ 


1041.293 318(239) 


0.203 559 482 1(89) 


4.483 28(93) 


-0.135(80) 




02 u 


1047.042 048(11) 


0.203 480 485(12) 


4.419 64(63) 


-0.712(73) 




Ol'l 


1372.459 242(136) 


0.202 657 042(22) 


4.542 71(63) 


[0.1] b 




03 3 


[1562.611 159] 


0.203 905 589 7(148) 


4.550 22(80) 


[-0.3] 




03'0 


1573.366 413(158) 


0.203 762 735 1(132) 


4.453 60(113) 


-0.801(206) 




00°2 


1710.976 247(76) 


0.201 635 352(32) 


4.533 64(111) 


0.396(102) 




02 2 1 


1886.947 787(135) 


0.203 048 230 7(117) 


4.639 14(113) 


[0.1] 




02°1 


1892.230 557(91) 


0.202 953 496 6(220) 


4.554 09(123) 


[-0.3] 




10°0 


2062.200 841(121) 


0.201 641 530 0(477) 


4.409 80(260) 


4.551(470) 


-0.1155(267) 


04"0 


[2084.378 374] 


0.204 248 031 9(608) 


4.592 0(296) 


[-0.3] 




04*0 


2099.524 648(251) 


0.204 051 914 6(232) 


4.494 89(828) 


[-0.6] 




WO 


2104.827 673(87) 


0.203 968 086 4(171) 


4.306 15(970) 


-32.73(814) 




01'2 


2218.028 446(175) 


0.202 091 366 5(527) 


4.682 71(123) 


[0.6] 




03 3 1 


[2402.340 630] 


0.203 428 111(148) 


4.734 93(580) 


[0.1] 




03' 1 


2412.122 352(193) 


0.203 259 077 6(203) 


4399 03(43) 


[-0.3] 




00°3 


2555.991 217(126) 


0.201 006 219(100) 


4.647 60(514) 


1.61(70) 




ll'O 


2575.307 586(181) 


0.202 015 427(45) 


4.473 84(210) 


2.86(28) 




05 5 


[2606.596 100] 


[0.204 589 3] 


[4.55] 


[-0.3] 




05 3 


2625.607 054(200) 


0.204 345 249(246) 


4.583 7(402) 


[-0.7] 




05'0 


2635.589 700(206) 


0.204 198 969(125) 


4.387 40(643) 


[-0.9] 




02 2 2 


[2726.564 984] 


[0.202 523 2] 


[4.804 7] 


[0.1] 




02°2 


2731.399 122(245) 


0.202 414 701(227) 


4.707(109) 


[-0.3] 




10°1 


2918.104 865(255) 


0.201 102 979(109) 


5.071 51(681) 


54.42(147) 


-1 .988(95) 


04 4 1 


[2918.572 190] 


[0.203 799 95] 


[4.34] 


[-0.3] 




04H 


293Z216 820(310) 


0.203 569 074(81) 


4.668 2(13) 


[-0.3] 




04"1 


2937.146 843(207) 


0.203 436 472(357) 


3.943 2(255) 


-52.44(634) 


1.71(50) 


01'3 


3057.093 032(518) 


0.201 510 88(145) 


4.824(81) 


[2.2] 




V?0 


[3088.908 46] 


0202 382 585(146) 


4.531 9(48) 


[2.2] 




12°0 


3095.554 42(9) 


0.202 311 240(147) 


4.500 9(56) 


[2.2] 




03 3 2 


[3236.415 907] 


[0.202 936 096] 


[4.776 6] 


[-0.2416] 




03 '2 


3245.260 572(251) 


0.202 744 75(42) 


4.771 0(128) 


[-0.2416] 




00"4 


3393.969 128(594) 


0.200 363 45(173) 


4.792(96) 


[2.5] 




11H 


3424.139 675(235) 


0.201 515 951(235) 


4.757 4(105) 


8.78(125) 




02 2 3 


[3560.106 489] 


[0.201 981 318] 


[4.983 28] 


[2.397] 




02"3 


3564.479 808(987) 


0.201 863 71(335) 


4.989(220) 


[2.397] 




13 3 


[3603.006 802] 


[0.202 745 521] 


[4.41] 


[-1.948] 




13 ! 


3615.345 30(20) 


0.202 617 950(637) 


4.569(66) 


6.7(189) 




04°2 C 


3762.825 61(83) 


0.202 601 23(931) 


15.06(323) 


-5459.(4814) 


160.(3450) 


10°2 d 


3768.497 40(22) 


0.200 840 86(229) 


-14.315(842) 


-3551.(1249) 


9138.(864) 


12 2 1 


3931.301 568(234) 


0.201 923 654(440) 


4.782 5(153) 


[6.00] 




12°1 


3937.427 356(330) 


0.201 835 302(539) 


4.745 3(180) 


[6.00] 




11'2 


4266.325 10(53) 


0.201 069 31(120) 


6.539 9(551) 


[240.0] 





* The uncertainty in the last digits (twice the estimated standard error) is given in parentheses. 
b The values enclosed in square brackets were fixed during the analysis. 

c Additional terms needed in the analysis were: M = - 1.03(118) x 10 -17 and N =4.26(153) x 10" 21 . See Rcf. [5.132] for discussion of 
analysis beyond / = 50. 
d Additional terms needed in the analysis were: M = -3.21(28) x 10 -17 and N =3.34(33) x 10~ 21 . 



430 



Volume 97, Number 4, July-August 1992 

Journal of Research of the National Institute of Standards and Technology 

Table 6. /-type resonance constants (in cm" 1 ) for 16 12 C 32 S 



Vib. state 
U1U2U3 



010 
020 
01 1 
030 
021 
040 
01 2 

03 1 
1 10 
050 
022 

04 1 
013 
120 
032 
1 1 1 
023 
1 30 
121 
1 12 



9v°Xl0 4 



2.121 938 68(53)' 
2.086 287(47) 
2.285 201(291) 
2.064 232(39) 
2.222 23(145) 

2.018 56(408) 
2.447 445(564) 
2.183 365(144) 
2.155 327(274) 

2.019 577(583) 
2.365(100) 
2.127 48(95) 
2.596 1(162) 
2.130 1(102) 
2.302 77(231) 
2.373 76(132) 

[2.42] 

2.100 35(154) 
[2.255 63] 

2.709 0(144) 



q vJ x 10" 



1.424 13(102) b 

0.659(35) 

3.593(114) 

0.252(22) 

1.898(64) 
- 12.49(345) 

6.146(177) 

1.475 6(245) 

3.764(70) 
-1.318(375) 

2.76(137) 
- 1.168(169) 

3.2(120) 

6.6(28) 

2.93(90) 

8.57(31) 

[0.6587]° 

2.06(87) 

[3.7] 
47.6(85) 



" The uncertainty in the last digits (twice the estimated standard error) is given in parentheses. 
b Also included in the fit was a higher order term q y jj = 0.574(44) x 10" IS . 
c The values enclosed in square brackets were fixed during the analysis. 



3.2.2 N2O The microwave data on N2O are 
not so extensive as for OCS but a great many mea- 
surements are still available. The early work of 
Pearson et al. [5.161] and of Lafferty and Lide 
[5.156] were very useful as were other early mea- 
surements given in the review by Lovas [5.5], 
namely the data given in Refs. [5.138, 5.139, 5.146, 
5.148, 5.159, 5.162]. The v\ rotational transitions 
given by Bogey [5.186] and the high-/ transitions 
given by Andreev et al. [5.194], by Burenin et al. 
[5.179], and by Vanek et al. [5.242] were of particu- 
lar value for better determining the centrifugal dis- 
tortion contribution to the line positions. 

All of the heterodyne measurements involving 
N2O have come from two laboratories, NRC in 
Canada, and NIST in the United States. The only 
saturated absorption measurements were those of 
Whitford et al. [5.183] on the laser transitions, 
10°0-00°1, near 930 cm -1 . The other infrared het- 
erodyne measurements were made by Wells and 
co-workers in a series of papers, Refs. [5.219, 
5.221, 5.224, 5.230, 5.231, 5.241, 5.243]. 

There have been a great many measurements on 
infrared bands of N 2 using either grating instru- 
ments or, more recently, FTS instruments. Some of 
the more important measurements which were 
used in the least squares refinement of the con- 



stants, but did not contribute to the band centers 
for the recommended calibration lines, were given 
in Refs. [5.191-5.193, 5.213, 5.215, 5.220, 5.225, 
5.229]. 

All these data were fit in the same way as was 
done for OCS. The constants given by the least- 
squares fit are given in Tables 7 and 8. 

For N2O the Fermi resonance is much more 
important than for OCS; nevertheless, the Fermi 
resonance was ignored in the fits and only the 
/-type resonance was included in the analysis. For 
some levels the Fermi resonance causes an effec- 
tive centrifugal distortion quite different from that 
of the ground state. 

In order to show the position of some of the hot 
band lines in the spectra, it was necessary to use 
constants for some levels not included in Tables 7 
and 8. The wavenumbers for those lines were 
calculated by taking the constants for the upper 
state reported in Refs. [5.191, 5.213, 5.215, 5.229]. 

For the less abundant isotopic species, the mi- 
crowave data given in Refs. [5.194, 5.227] were 
used in preparing these tables. One paper has re- 
ported heterodyne frequency measurements in the 
infrared for 15 N u N lfi O and 14 N 15 N lfi O [5.231]. For 
the most part the data for the rarer isotopic species 
were taken from Refs. [1.7, 5.173, 5.184, 5.192, 
5.193, 5.213, 5.215, 5.225, 5.229]. 
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Table 7. 


Rovibrational constants (in cm" 


"') used for "N'WO 








Vib. state 


M) 


B, 


DsrXlO 7 


/f v x 10 14 


LyXW 17 


WO 


0.0 


0.419 011 006(15)* 


1.760 91(19) 


-1.66(21) 




Ol'O 


588.767 741(163) 


0.419 573 590(23) 


1.788 70(33) 


-0.82(46) 




02°0 


1168.132 418(198) 


0.419 919 855(52) 


1.871 21(85) 


-8.18(368) 




02*0 


1177.744 555(85) 


0.420 124 817(41) 


1.816 72(78) 


5.25(504) 


-2.14(69) 


00°1 


1284.903 289(124) 


0.417 255 066(20) 


1.725 67(24) 


11.30(42) 


0.444(27) 


03'0 


1749.064 972(166) 


0.420 331 191(67) 


1.911 42(80) 


-11.89(205) 




03*0 


1766.911 896(161) 


0.420 664 523(346) 


1.852 1(64) 


17.4(405) 




Ol 1 ! 


1880.265 695(150) 


0.417 918 446(50) 


1.733 61(53) 


14.15(126) 


0.288(102) 


10°0 


2223.756 693(124) 


0.415 559 512(18) 


1.754 67(20) 


-1.359(218) 




04°0 


2322.572 934(211) 


0.420 618 036(216) 


1.943 02(155) 


-384.4(164) 




04^ 


2331.121 460(124) 


0.420 768 166(253) 


1.990 06(230) 


321.2(158) 




04^ 


[2356.251 397f 


[0.421 193 718] 


[1.90] 


[0.0] 




02 u l 


2461.996 447(242) 


0.418 147 317(318) 


1.892 76(175) 


[0.0] 




02 2 1 


2474.798 428(324) 


0.418 530 238(655) 


1.750 31(257) 


[0.0] 




00°2 


2563.339 334(169) 


0.415 605 588(268) 


1.639 20(152) 


64.40(224) 




ll'O 


2798.292 466(200) 


0.416 159 111(81) 


1.782 02(46) 


-0.47(58) 




03'1 


3046.212 560(931) 


0.418 567 39(164) 


1.904 70(539) 


[0.0] 




03 3 1 


[3068.720 525] 


[0.419 107] 


[1.818] 


[0.0] 




01'2 


3165.853 959(214) 


0.416 382 245(410) 


1.632 92(294) 


48.70(571) 




13'0 


3931.248 302(341) 


0.416 994 90(118) 


1.909 53(549) 


[0.0] 




13 3 


3948.285 330(370) 


0.417 327 81(195) 


1.835 5(123) 


[0.0] 





1 The uncertainty in the last digits (twice the estimated standard error) is given in parentheses. 
b The values enclosed in square brackets were fixed during the analysis. 



Table 8. /-type resonance constants (in cm" 1 ) for 14 N 14 N 16 



Vib. state 



^"xlO 4 



g VJ x 10 9 



10 
020 
030 

01 1 
040 
021 
1 10 
03 1 
12 
130 



7.920 055 2(83)* 
7.607 30(98) 
7.472 506(257) 
9.083 842(493) 
7.481 95(567) 
8.206 5(118) 
7.771 651(471) 
8.086 4(103) 
10.722 70(197) 
7.326 0(111) 



1.002 0(155) 
2.766(108)" 
2.889(49) c 
-2.877(32) d 
12.379(424) 
1.405(305) 
1.199 1(185) 
2.292(416) 
- 10.713(76) 
3.062(567) 



"The uncertainty in the last digits (twice the estimated standard error) is given in parentheses. 
b Also included in the fit was a higher order term q VJJ = 0.123(315) x 10" 13 . 
'Also included in the fit was a higher order term q w ,,= -0.224(102)xl0"' 3 . 
d Also included in the fit was a higher order term q vlj = 1.109(32)- 10"' 3 . 



3.2.3 CS 2 Since carbon disulfide (CS 2 ) is a 
symmetric linear molecule, it is nonpolar and has 
no microwave spectrum. There are, however, sev- 
eral high resolution infrared studies of its spectrum 
in the 1450 to 1550 cm -1 region [5.383, 5.386, 
5.389] in addition to the heterodyne measurements 
made by Wells et al. [5.390]. A number of other 
measurements have been made on CS2 so that the 
ground state constants, B a and Do, are quite well 
determined for both the 12 C 32 S Z and 13 C 32 S 2 iso- 
topic species. 



The analysis of the V3 band of CS2 is uncompli- 
cated by either Fermi resonance or /-type reso- 
nance. The analysis used Eqs. (3.1)-(3.3) as 
described in detail by Wells et al. [5.390]. The 
recommended calibration frequencies are based on 
the constants given in their paper. Only the 
00°1-00°0 transitions of 12 C 32 S 2 and "C^Sa should 
be used for calibration. The other line wavenum- 
bers given in the tables are for identification pur- 
poses and to show how close some weaker lines 
may be to the calibration lines. The wavenumbers 
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of the other lines were calculated from the con- 
stants given by Winther et al. [5.389] and may be in 
error by as much as 0.01 cm -1 . In the spectral maps 
some lines may not be shown in the region below 
1505 cm -1 because they arise from transitions not 
included in the data base. 

3.2.4 CO The carbon monoxide (CO) 
wavenumbers for the calibration lines given in the 
atlas were calculated from the constants given in 
Table 9 and based on Eqs. (3.2) and (3.3). Only the 
wavenumbers given for the 12 C 16 molecule were 
determined adequately by frequency measure- 
ments and so they are the only wavenumbers that 
should be used for frequency (wavenumber) cali- 
bration. There are a few frequency measurements 
for the other isotopomers, but not enough to 
provide good calibration. 

For the 12 C 16 molecule the ground state con- 
stants are primarily based on the sub-millimeter 
wave measurements given by Gordy and Cowan 
[5.250], Rosenblum et al. [5.253], and Helminger et 
al. [5.266] and on the far-infrared heterodyne mea- 
surements of Nolt et al. [5.309] and Varberg and 
Evenson [5.320]. Also included in the fit were mi- 
crowave measurements of the / = l<-0 transitions 
in the first two vibrationally excited states as 
reported by Dixon [5.288]. Aside from these two 
measurements the upper state constants are based 
primarily on the heterodyne measurements of the 
1-0 band given by Schneider et al. [5.314] and by 
Maki et al. [5.316] and the 2-0 band given by 
Pollock et al. [5.304]. Also used was one sub- 
Doppler measurement of the 1-0 band communi- 
cated to us by Urban [5.321]. 

Other data included in the analysis to determine 
the best constants were some heterodyne laser 
measurements given by Schneider et al. [5.313] and 
some FTS measurements from Guelachvili et al. 
[5.289, 5.301, 5.307] and Brown and Toth [5.306]. 
The FTS measurements were used to help deter- 
mine the best centrifugal distortion constants for 
the v = l, 2, and 3 states. Only heterodyne fre- 
quency measurements were used to determine the 
vibrational frequencies. 

The analysis of the CO data was carried out in 
the same way as described by Maki et al. [5.316] 
except that new data have been included [5.320, 
5.321]. In order to avoid the possibility of problems 
with the potential function model, the lowest order 
constants, Yio, Y20, Yin, Yn, Y21, Y02, Y12, and Y03 were 
fit to data for only the v = 0, 1, and 2 states. The 
other constants were constrained to values given by 
earlier fits which included data for higher vibra- 



tional states, as described by Schneider et al. 
[5.314]. 

The uncertainties given for the calibration 
wavenumbers are based on the variance-covariance 
matrix given by the least -squares fit with the higher 
order constants constrained. 

The wavenumbers for the other isotopic species 
of CO were calculated from constants given by 
Guelachvili et al. [5.301], but corrected to agree 
with the offset observed in the wavenumbers for 
the 12 C 1(i O species. The wavenumbers for the rarer 
isotopic species are given to help in correctly iden- 
tifying the calibration lines. 



Table 9. Constants used to calculate the u = l«-0, t» = 2«-0, 
and u = 2<-l transition wavenumbers for ,2 C ,6 



Constant" 



Wavenumber b 
(cm" 1 ) 



Y*. 

Y M 

Y40XIO 4 

YjoXlO 6 

Yo, 

Y„ 

Y21XIO 6 

Y 3 , x 10 7 

YkxIO 5 

Y^xlO 8 

Y22X10" 

Y03XIO 11 
Y, 3 xlO' 2 
YhXIO 16 
Y, 4 xl0' 8 
YojxIO 22 



2169.812 
- 13.287 

0.010 

0.740 
-0.137 

1.931 
-0.017 

0.486 

0.333 
-0.612 

0.100 
-0.177 

0.589 
-0.145 
-0.360 
-0.684 
-0.471 



615(26) c 

812 0(87) 

383 46(980) 

03(1300) 

37(5913) 

280 858 2(555) 

504 036 7(1302) 

50(4041) 

87(2754) 

159 11(230) 

669(3765) 

14(2551) 

265(709) 

467(300) 

976(20) 

5(47) 

36(167) 



* Dunham coefficients defined by Eq. (3.2). 

b To convert to frequency units multiply by 29 979.2458 MHz/ 

cm -1 . 

"The uncertainty (twice the estimated standard error) in the 

last digits is given in parentheses. Excess digits are given to 

avoid round-off errors. 



3.2.5 NO The wavenumbers of the line posi- 
tions used to produce the NO atlas were calculated 
using the constants given by Hinz et al. [5.371] for 
the 14 N 16 species. The constants given by Amiot 
and Guelachvili [5.355] were used for 15 N 16 and 
those given by Amiot et al. [5.350] were used for 
itjsjisQ Tho Se were the only transitions strong 
enough to show on the spectral plots. The uncer- 
tainties given in the tables are only estimates based 
on the accuracy of other heterodyne measurements 
and an estimate of the accuracy of the Hamiltonian 
used to fit the data. 
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The tables only give the wavenumber values for 
lines of 14 N 16 0. Because only one paper [5.371] re- 
ports infrared frequency measurements for NO, 
the NO lines are not recommended for calibration 
in those regions where either OCS or N 2 transi- 
tions are available for calibration. For many spec- 
trometers most of the NO transitions are 
unresolved doublets. At low pressures the doublets 
will have equal intensities and widths so that no 
appreciable errors will be incurred by using the av- 
erage frequency of the doublet. A few of the low-/ 
transitions, especially for the Q -branch lines, will 
have additional structure due to the interaction of 
the nuclear electric quadrupole moment of the ni- 
trogen atom with the surrounding charge distribu- 
tion. Such small splittings will be no larger than a 
few megahertz and the average line positions given 
in the tables will not be affected. 

33 Intensity Calculations 

This section is intended to show how the transi- 
tion intensities were calculated and how the inten- 
sity was defined. This will enable users to 
determine the intensity to be expected for condi- 
tions other than those used to prepare the atlas. In 
order to estimate the appearance of the spectrum 
under conditions of pressure broadening and spec- 
trometer resolution different from that used to 
produce the atlas figures, it is also necessary to 
consider the line shapes and the effect of finite slit 
functions. Such effects are discussed in Sec. 4. 

The integrated intensity of an individual line 
representing a single rovibrational transition is in- 
dependent of the line shape. For this atlas we take 
the integrated line intensity, 5, to mean 

S=l k(v)dv = (l!pl)\ ln(h/I)„dp, (3.9) 

J — 00 J — 00 

where k(v) is the absorption coefficient at fre- 
quency v, p is the partial pressure of the gas, / is 
the length of the absorption path, / is the intensity 
of radiation without absorption from the line in 
question, and / is the intensity of radiation after 
absorption by the line. Note that Eq. (3.9) is only 
concerned with the absorption due to a particular 
transition. 

The integrated intensities, as given in this atlas, 
were calculated by using the equation 

5 = exp( -E"lkT)[\ - exp( - W0.69504T)] 

(WG.Gr) X vC\Rffl£? S 2 S£ , (3.10) 



where C is a proportionality constant that includes 
the factors 8ir 3 /3hc and other factors, such as the 
Loschmidt constant (2.686 763 x 10 25 molecules/ 
m 3 ), required to give 5 in appropriate units. If S is 
in units of cm/molecule at STP, C is 4.162 
38 x 10 -19 cm 2 D -2 /molecule. If 5 is in units of cm -2 
atm" 1 at temperature T, C would be 3054.7262/7 
cm -1 D" 2 atm -1 . In Eq. (3.10) \R$I%\ is the vibra- 
tional transition moment or dipole derivative in de- 
bye units (ID = 3.335 64xl0 _30 Cm), v is the 
transition wavenumber in units of cm -1 , T is the 
temperature in kelvin, N-, is the concentration of 
the isotopic species under consideration, Q v and Qr 
are the vibrational and rotational partition func- 
tions respectively, and S v and 5r are vibrational and 
rotational strength factors that should be included 
in the intensity. Some workers prefer to include 5 V 
in the transition moment but we prefer to express it 
separately so the transition moment can be seen to 
be nearly the same for the ground state transitions 
and the accompanying hot bands. 5r is also called 
the direction cosine matrix element. The intensities 
given in these tables were calculated for a tempera- 
ture of 296 K. 

In Eq. (3.10) the term Ni/QvQr compensates for 
the fact that the pressure used in Eq. (3.9) is the 
total pressure of the gas being measured, including 
all isotopic species. That is to say, the pressure 
does not take into account the isotopic concentra- 
tion or the number of molecules in different states. 
The values of N-, were calculated from the isotopic 
abundances given in Table 10 and taken from Refs. 
[3.5, 3.6]. The vibrational partition function was 
calculated by summing the Boltzmann population 
of the vibrational energy levels. The vibrational 
partition functions are given in Table 11. 

The rotational partition functions were calcu- 
lated from the equations given by McDowell [3.7]. 
A different rotational partition function was calcu- 
lated for each vibrational energy level and for each 
isotopomer. 

Some workers like to use S°=S/Ni for 5 because 
5° seems to be a more appropriate molecular 
property. On the other hand, S is more useful for 
analytical purposes, such as the determination of 
the amount of CO in the atmosphere. Actually, it is 
\R\ or \R\ 2 , rather than 5, that is the true molecu- 
lar property in Eq. (3.10), and, as suggested by 
Toth [5.218], all intensity measurements should re- 
port the value of \R\, or \R\ 2 . Unfortunately, some 
authors have left out the Ni term and report values 
of \R\ that are not true molecular properties relat- 
able to the electron distribution in the molecule. 
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In the infancy of infrared spectroscopy the in- 
strumentation was unable to resolve individual 
rovibrational transitions, so many early papers 
measured the intensity of entire vibrational bands. 
Those papers measured band intensities by using a 
modification of Eq. (3.9) in which the integration 
was over the entire band rather than over a single 
line. The fine points of isotopic concentration and 
vibrational hot bands were ignored. Even today in- 
tegrated band intensities are often measured for 
heavier molecules for which the density of lines is 
very high. 

From those early measurements the term band 
intensity came to mean the intensity of all the lines 
in a band including all isotopes present in a normal 
sample and all hot bands. As a first approximation 
such a band intensity can be calculated from Eq. 
(3.10) if the terms S?S 2 R , and NJQ V Q R are all omit- 
ted and v is set equal to the center of the band. 
This band intensity will generally be within a few 
percent of the intensity obtained by adding all the 
line intensities for the band. To be faithful to the 
original meaning, the true band intensity should be 
the sum of the intensities of all the lines in the 
band, including all isotopes and all hot bands. 

Table 10. Percent isotopic abundances 



12 C 
»C 

14 N 
I5 N 
16Q 

"O 
18Q 

32 s 
33 s 

34 S 

36 S 



98.90 
1.10 

99.634 
0.366 

99.762 
0.038 
0.200 

95.02 
0.75 
4.21 
0.02 



Table 11. Vibrational partition function and natural abundance 
for the various isotopic species found in the calibration atlas 



Molecule 


g v (296K) Ni 


Molecule Q v (296 K) N { 


»o 12 c? 2 s 


1.199 


0.937 5 


12 c lfi o 


1.000 


0.986 6 


i6 o 12 c? 4 s 


1.203 


0.041 63 


13 C 16 


1.000 


0.010 97 


i6 o 13 c 32 s 


1.216 


0.010 50 


12 C 18 


1.000 


0.001 98 


16 12 C 33 S 


1.201 


0.007 40 


12 c 17 o 


1.000 


0.000 38 


i8 o 12 c 32 s 


1.207 


0.001 92 


13 C 18 


1.000 


0.000 022 


14 N 14 N 1S 


1.1273 


0.990 3 


13 c 17 o 


1.000 


0.000 004 


14 N 15 N W 


1.1364 


0.003 64 


12 c 32 s 2 


1.424 


0.892 95 


15 N 14 N lfi O 


1.1298 


0.003 64 


13 C? 2 S 2 


1.452 


0.009 93 


14 N 14 N 1« 


1.1308 


0.001 99 


12 c 34 s 32 s 


1.430 


0.079 13 


14 N 16 


1.000 


0.994 


12 C 33 S 32 S 


1.428 


0.014 10 


"N lfi O 


1.000 


0.003 7 








14 N 18 


1.000 


0.002 








14 N 17 


1.000 


0.000 38 









3.3.1 Calculation of the Strength Factors 

The vibrational strength factors were slightly dif- 
ferent for even and odd values of Avi. For simplicity 
the strength factor was broken into two factors such 
that 



S v — 5i352 



(3.11) 



with 

5, 2 3 = (u, + 4u 1 )! (vi + Aviy/WvilAvilAvil), (3.12) 
and for Av2 even (or zero) 

Si = [1/2(^2 + / + Av2)]\ [1/2(^2 - / + Avi)][ 

/{[V2(v2+l)]\ [1/2(iW)]!x{[1/2(4v2)]!} 2 }, (3.13) 
while for Av2 odd 
Si = [1/2(U2 + / + Avj - 1)]! [1/2(U2 -/ + Av2 - 1)]! 
(V2+IAI + Av2 + 1)/{[1/2(u2+/)]![1/2(u2-/)]! 



{[l/2(^l^-l)]!} 2 (4^ + 1)}. 



(3.14) 



In both cases 5 2 =0 if |A/| > 1. In Eqs. (3.12)-(3.14) 
the smaller of v' and v" is used for v, Av = \v' — v"\, 
I =/", and A/ =/'-/". For the present calculations, 
52 can be taken as the positive square-root of Eqs. 
(3.13) and (3.14). 

Equation (3.12) was derived from the properties 
of harmonic oscillator wave functions such as can 
be found, among other places, in Appendix III of 
Ref. [3.8]. Equations (3.13) and (3.14) were 
derived from the properties of the two-dimensional 
harmonic oscillator wave functions given by Moffitt 
and Liehr [3.9]. 

Note that S 2 is normalized so that transitions 
from the ground state always have S 2 = l. For cer- 
tain hot bands, such as 2^i - v\, S 2 = 2, while for the 
hot band 7>v\. — v\., 5 2 = 3. Other authors sometimes 
include 5 V in the transition moment, \R\, in which 
case the transition moment for certain hot bands 
will be very different from the transition moment 
for the ground state transitions. 

For transitions for which t>2 = (and / = 0) the 
rotational strength factors are given by 



and 



S R = \m\ for AJ = ±1 
Sk = 0foiAJ = 0, 
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where m has the usual meaning of — /" for 
AJ=-1 and/" + l for AJ= + 1. 

If 1*2^0, the /-type resonance energy matrix was 
used and the intensity was obtained by multiplying 
each term of the eigenvector, for the appropriate 
eigenvalue, by the appropriate intensity factor 
given by S v Sr where Sr was determined from Table 
2.1 of Gordy and Cook [3.11], or Table 4-4 of 
Townes and Schawlow [3.12]. Di Lauro and Mills 
[3.13] describe a similar procedure for determining 
intensities of transitions in Coriolis coupled levels 
of a symmetric rotor. For the specific cases of /- 
type resonance, this procedure was described by 
Maki et al. [3.14] although they dealt with the sym- 
metry factored matrix whereas the present calcula- 
tions used the unfactored matrix, such as Eq. (3.7). 

To understand this intensity calculation let us 
consider a transition from an unperturbed lower 
state (ui,i»2,/,U3,/) to an upper state that is in- 
volved in /-type resonance, (vi',vi*,l ',Vi',J') = (i»i + Avi, 
vi + Avi,l + Al,Vi+ AviJ '). The upper state energy is 
given by a particular eigenvalue, E\. The eigenvec- 
tor for this eigenvalue gives the mixing coefficients 
an, «i2, ai3, etc. that measure the contribution of 
each unperturbed state to the perturbed state. The 
transition intensity is given by 

| < i»i + Av\,vz + AviJ. + A1,V3 + Av 3 ,J + AJ\/jl\ 

vi,v 2 ,l,v 3 ,J>\ 2 = 

|fln < vi + Avi,vi + Avi,l ' = vi + Avi,V3 + Av3,J 

+ Aj\(ju\vi,V 2 ,l,V3j> 

+ an<Vi + Avi,V2+ Av2,l' 

= i»2 + Av2 — 2,i»3 + AviJ + AJ | \x, | i»i,i»2,/ ,i»3,/ > 



If one assumes that the transition moment, or 
dipole derivative, is the same for all values of /, 
then Eq. (3.15) can be rewritten 

\R fSlSl = \R \ 2 Sh {auS 2 S R (v2 + Av2,l 

= V 2 +Av 2 ,J'<r-V2j",J")+ai2S2S R (v2 + Av2j 

= V2 + Av2-ZJ'<r-V2j",J") +a 13 . . .} 2 . (3.16) 



Each term within the curly brackets on the right 
side of Eq. (3.16) has a different value for S2 and 
5r depending on the unperturbed transition to 
which they apply. In many cases they will be zero 
because they apply to |4/|> 1. In Eqs. (3.12)-(3.14) 
we have already given the formulas for the values 
of Sh and Si 

The only nonzero values of Sr are given below. 
For Al = 0, the AJ = transitions have 



5 R = [(2/ + l)/ 2 //(/ + l)] 1/2 
and the AJ = ± 1 transitions have 



(3.17) 



5 R = [(|m| 2 -/ 2 )/|m|] w . (3.18) 

For Al = ±\, the expressions for AJ = were 

5 R = l/2[(2/ + 1)(/ + /)(/ -/ + 1)//(/ + 1)] 1/2 , (3.19) 

where / is the larger of / ' and /". For Al = + 1 and 
AJ= + 1, 

5 R = - 1/2[(/ +/)(/ +/ + 1)//] 1/2 ; (3.20) 



+ ai3<t>i + Avi,V2 + Av 2 ,l' 



= V2 + AV2 — 4,l>3 + AV3J + AJ\ fji\vi,V2,l ,V3,J > 



••I 2 - 



(3.15) 



All the upper state quantum numbers on the left 
side of Eq. (3.15) are the same as those on the right 
side except the /' values. The /-type resonance 
mixes levels that differ only in the value of/. 



for Al = + 1 and AJ = - 1, 

Sr = 1/2[(-/+/)(-/+/ + 1)//] 1/2 ; (3.21) 

for Al = - 1 and AJ = + 1, 

5 R = l/2[(/-/)(/-/ + l)//] 1/2 ; (3.22) 
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and for Al - - 1 and AJ=—\, 



Sk=~ l/2[( -/ -/)(-/-/ + 1)//] 1/2 . (3.23) 

In Eqs. (3.19)-(3.23) / is the larger of/' and J" but 
/ is always /". 

When both the upper and lower states are in- 
volved in /-type resonance, Eqs. (3.15) and (3.16) 
must be modified to include the eigenvectors for 
both the upper and lower states, otherwise the 
treatment is the same. 

3.3.2 Herman-Wallis Terms Sometimes it is 
necessary to divide \RirPr\ 2 into two terms such that 



|Di)T/'|2_|pi)V|2 E-/' 

\Kv"rr\ — \tiv"r\ ?r , 



(3.24) 



where Fr is similar to the Herman-Wallis term 
[3.15, 3.16]. There are several forms that have been 
used for the Fr term. For the overtone of CO the 
form used was 



F r = l + Cim+C 2 m 2 . 



(3.25) 



Toth [5.218], in his extensive intensity measure- 
ments for N 2 0, has given the values for the 
Herman-Wallis constants for several bands. Toth 
has used several formulations such as Eq. (3.25) 
above and also 

Ff = [l+aJ(J + l)+bJ 2 (J + l) 2 ] (3.26) 



F/ = [l+flim+fl 2 /'(7' + l)] 2 (3.27) 



F/ = [l + fim+£>mT. 



(3.28) 



In the case of OCS, Dang-Nhu and Guelachvili 
[5.104] used the form of Eq. (3.25) for the ll'O- 
00°0 band. The same constant was also used to cal- 
culate the intensities of the hot bands for which 
At>i = 1 and Au 2 = 1. The same form and, coinciden- 
tally, almost the same value for the Herman-Wallis 
term was used for the vi band system. As recom- 
mended by Maki et al. [5.132], we used 



Ff = [l+a 3 /'(/' + !)] 



(3.29) 



for the hot bands 10°2-00°1 and 04°2-00°l in order 
to allow for the resonance interaction of the upper 



states. For the other transitions, there seemed to 
be no need for using a Herman-Wallis term in the 
intensity calculations. 

For the (B'O-Ol'O band of OCS, Depan- 
nemaecker and Lemaire [5.118] found that the first 
term, G, in Eq. (3.25) was needed although no 
Herman-Wallis term was needed for the 02°0-00°0 
band. On the other hand, the more extensive mea- 
surements of Blanquet et al. [5.134] showed that a 
Herman-Wallis term would improve the intensity 
fit for 02°0-00°0. 

3.4 Data Sources for Intensity Calculations 

3.4.1 OCS Only a small number of intensity 
measurements have been published for OCS. Most 
of them are measurements of integrated band in- 
tensities including transitions from the ground 
state as well as from other low energy states popu- 
lated at room temperature, hot bands. The most 
thorough studies of band intensities are those of 
Foord and Whiffen [5.49] and Kagann [5.89] where 
references to earlier work may be found. Kagann's 
values were consistently smaller by 11 to 15 percent 
except for the strongest band near 2062 cm -1 for 
which Kagann found an intensity that was 18 per- 
cent larger. 

There have been only nine papers reporting in- 
tensity measurements of individual rovibrational 
transitions [5.73, 5.74, 5.83, 5.85, 5.103, 5.104, 
5.118, 5.127, 5.134] and the first three of those pa- 
pers are considerably less thorough than the oth- 
ers. In addition, there are two papers that give 
information on the relative intensities of well re- 
solved rovibrational transitions without making ab- 
solute intensity measurements [5.101, 5.132]. 

Several papers [5.49, 5.97, 5.100, 5.106, 5.112] 
have used intensity and/or dipole moment mea- 
surements to derive a dipole moment function for 
OCS. At the present time the dipole function of 
OCS can be used to predict the dipole moment of 
OCS in different vibrational states, but it is not 
very useful for predicting transition intensities. The 
accuracy of the dipole moment function for pre- 
dicting dipole moments is due largely to the fact 
that it is based on very accurate laser-Stark mea- 
surements involving many vibrational states. 

Table 12 summarizes the intensity measurements 
reported in the literature. Table 12 also gives the 
transition moments used in the calculations for this 
atlas and the integrated intensity of each band as 
obtained by actually summing all of the transitions 
within each band, including all hot bands and all 
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isotopic species. In a few cases, where indicated by 
a footnote, the integrated intensity given in the lit- 
erature and reported in Table 12 does not include 
the hot band intensity. In those cases Table 12 
gives the intensity calculated without including the 
contribution from hot bands in order to more easily 
compare the intensity used in this atlas with what 
was reported in the literature. 

Only a few bands of OCS show any need for 
Herman-Wallis terms. As is usually the case, the 
two perpendicular bands, vi and vi + i>2, required a 
small Herman-Wallis term, d = 0.0045 and 0.00463 
for the two bands, respectively. The Herman-Wal- 



lis term is defined by Eqs. (3.24) and (3.25). For n 
it was estimated by us from inspection of the exper- 
imental spectrum and is not a very accurate value. 
The Herman-Wallis term for v\ + vi was measured 
by Dang-Nfhu and Guelachvili [5.104]. Although 
the Herman-Wallis term is likely to be somewhat 
different for different isotopes and for hot bands, 
the same term was used for all bands involving the 
same quantum number changes unless indicated in 
the key that accompanies the tables. 

The only other bands that were given nonzero 
Herman-Wallis constants were some of the hot 
bands that go with 2^2 and the hot bands 04°2-00°l 



Table 12. Integrated intensities and transition moments for OCS 











Used in Atlas Calculation 




Band 


Frequency 


Measured 


Ref. 


Integrated 




Transition 




interval 


integrated 




intensity 




moment" 






intensity 




(cm _2 atm') (cm/molecule) 


IW 




(cm" 1 ) 


(cm- 2 atm _I ) c 
at 296 K 




at 296 K 


X10 1 " 


(D)« 


V2 


494-567 


11.9(8) c 
12.0(8) 


[5.28] 
[5.89] 


11.85 


4.78 


0.047 


Vi 


812-S90 


41.(4) 

35.4(16) 

35.5(23) 

36.3(10) 

36.3(1) 


[5.49] 

[5.85] 

[5.89] 

[5.103] 

[5.127] 


36.3 


14.64 


0.064 


2U2 


1000-1095 


14.9(15) 
13.3(9) 
12.9 r 
12.7 f 


[5.49] 
[5.89] 
[5.118] 
[5.134] 


13.8 f 


5.57 


0.0333 


2^ 


1650-1739 


7.6(16) 
6.71(47) 


[5.49] 
[5.89] 


6.9 


2.78 


0.0195 


2vi+ i>i 


1832-1934 


12.4(25) 
10.9(8) 


[5.49] 
[5.89] 


11.2 


4.52 


0.022 


"i 


1970-2091 


2520.(250) 


[5.49] 


2533. 1021.7 


0.345 






2980.(220) 


[5.89] 








4^ 


2080-2141 


19.5(49) 
9.6 f 


[5.49] 
[5.101] 


10.4 ' 


4.19 


0.024 


3^ 


2510-2574 




[5.132] 


035 


0.14 


0.00367 


»i + vi 


2550-2600 


0. 147(6) f 


[5.104] 


0.149' 


0.060 


0.00257 


2o2+2vj 


2693-2763 


0.52(8) 


[5.49] 


0.52 


0.21 


0.004 


V, + V3 


2862-2941 


37.6(40) 
33.6(22) 


[5.49] 
[5.89] 


35.5 


14.3 


0.034 


4j>2+l>) 


2910-2970 


0.32 f 


'[5.132] 


0.32 f 


0.13 


0.0036 


1>1+2j>2 


3065-3120 


3.01(30) 
2.25(16) 


[5.49] 
[5.89] 


2.43 


0.98 


0.008 



* See the key page for each band region for transition moments of hot bands that are not the same as for the ground state 

transitions. 

b See Eq. (3.24) for definition of \RX-l-\. 

c 1 cm -2 atm -1 at 296 K=4.0335x 10" 20 cm/molecule. 

<, lD = 3.336xlO- M Cm. 

e The uncertainties in the last digits are given in parentheses. 

'Hot band intensities were not included. 
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and 10°2-00°1 which involve resonance coupled up- 
per states. Depannemaecker and Lemaire [5.118] 
determined that C\- 0.0019 gives the best fit to 
their intensity data for (B'O-Gl'O transitions. We 
have used the same Herman- Wallis term for most 
of the other hot bands accompanying 2vi\ see the 
descriptive key accompanying the 1000 to 1095 
cm -1 tables. Maki et al. [5.132] have given effective 
constants a 3 = 0.00105 and a 3 = -0.00034, respec- 
tively, to the resonance coupled transitions. The 
latter constants are based on Eqs. (3.24) and 
(3.29). 

After the tables were prepared the new intensity 
measurements of Blanquet et al. [5.134] became 
available and it leads us to believe that the inten- 
sity values in the tables are too high by about 7 



percent because we assumed that the dipole 
derivative given by Depannemaecker and Lemaire 
took into account the isotopic abundance. Any re- 
vised intensity calculations for the 2vi band should 
also include the Herman- Wallis constant measured 
by Blanquet et al. [5.134]. 

3.4.2 N2O A great many intensity measure- 
ments have been made on the various bands of N2O 
given in this atlas. Not all the measurements are 
given in Table 13, but the most important or most 
recent measurements are given there. 

As was the case for OCS, the earlier measure- 
ments were of the intensity of the unresolved bands. 
Such measurements included all hot bands and all 
isotopic species present in a normal sample. One of 



Table 13. Integrated intensities and transition moments for N2O 











Used in Atlas Calculation 




Band 


Frequency 


Measured 


Ref. 


Integrated 




Transition 




interval 


integrated 




intensity 




moment" 






intensity 




(cm -2 atm 1 ) (cm/molecule) 


IW 




(em- 1 ) 


(cm -2 atm -1 )' 
at 296 K 




at 296 K 


X10 19 


(D) d 


V2 


520-660 


28.9(15) c - f 
30.4(40) ' 
31. 4(40) l 


[5.214] 
[5.151] 
[5.157] 


24.5 
29.2 f 


9.88 


0.0692 


Vl—Vs 


880-990 


0.0524 
0.055 


[5.209] 
[5.172] 


0.052 


0.021 


0.056 


V,-2V2 


990-1090 


-- 


-- 


0.01 


0.004 


0.018 


2^2 


1110-1225 


7.00(3) 
7.36(7) 
6.98 

8.27(43) < 
8.49(9) f 


[5.218] 
[5.216] 
[5.237] 
[5.214] 
[5.234] 


6.66 
8.30 f 


2.69 


0.0247 


V3 


1200-1340 


207.(1) 
225.(12) f 
226.(2) f 
229.8(45) 


[5.218] 
[5.214] 
[5.234] 
[5.216] 


206. 
232.5 f 


83.1 


0.1326 


Vl+Vi 


1835-1925 


0.52(3) 

0.58(10)' 

0.437(20) 


[5.212] 
[5.214] 
[5.185] 


0.52 
0.60 f 


0.21 


0.00548 


Vl 


2140-2269 


1207.(22) 

1421.(76) f 

1189.(30) 

1302.(50) 

1301.(54) f 


[5.226] 
[5.214] 
[5.187] 
[5.171] 
[5.198] 


1206. 
1358. f 


486.4 


0.244 


2xi+^3 


2410-2510 


6.79(4) 
7.35(38) f 


[5.216] 
[5.214] 


6.7 
8.2 f 


27 


0.017 


2vi 


2490-2605 


30.86(14) 
32.1(17) f 


[5.216] 
[5.214] 


32.1 
36.3 < 


12.9 


0.037 


V1 + V2 


2725-2840 


2.24(13) f 
2.11 


[5.214] 
[5.235] 


1.895 
2.26 f 


0.764 


0.0086 



" Unless indicated otherwise (see footnote f), the integrated intensities are given only for the vibrational transition from the 

ground state. 

b See Eq. (3.24) for definition of |R#'|. 

c 1 cm -2 atm -1 at 296 K=4.0335 x 10 -20 cm/molecule. 

d lD=3.336xl0- 3O Cm. 

c The uncertainties in the last digits are given in parentheses. 

r All hot bands and isotopes are included. 
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the more complete sets of measurements is that 
reported by Kagann [5.214]. Many earlier measure- 
ments are summarized in that paper. Kobayashi 
and Suzuki [5.228] have fit the available data on 
intensities to a dipole moment function from which 
they calculate transition moments and Herman- 
Wallis constants. More measurements are needed 
to evaluate the accuracy of the constants calculated 
from their dipole moment function. 

Because of the interest in N2O as an atmospheric 
gas, the more recent intensity measurements in- 
volve spectra with resolved rotational structure and 
individual rovibrational lines have been measured. 
In many places in Table 13 we give both the inten- 
sity of the transitions from the ground state and the 
intensity of all transitions with the same change of 
quantum numbers regardless of the lower state. 
That allows one to compare our intensities with 
those reported in the literature. 

For most of the transitions we have tried to use 
what we judge to be the best measurements for the 
intensity but in a few cases we have departed from 
the literature values in order to give a more realis- 
tic appearance to the spectrum. For example, we 
have used the transition moment given by Toth 
[5.218] for the vs band near 1270 cm -1 , but the 
transition moment given for the 2i^ band did not 
give a calculated spectrum that matched the ob- 
served spectrum where the two bands overlap. This 
mismatch caused us to use a slightly smaller value 
for the transition moment for 2i^ than that given by 
Toth. It is possible that either the Herman-Wallis 
constants or the intensity of the v$ band should 
have been changed. The more reasonable thing 
seemed to be to change the intensity of the weaker 
band. We have also recognized that Toth has in- 
cluded the vibrational factors in the dipole moment 
matrix elements that he reported, whereas Eq. 
(3.10) treats that as a separate term that multiplies 
the transition moment. 

The two bands 2i>j+ ^ and 2v$ also overlap and 
again there was some modification of the literature 
values for the intensities in order to get a good 
agreement between the calculated spectrum and 
that obtained experimentally. 

We know of no intensity measurements for the 
v\—2vi transitions between 990 and 1090 cm -1 ; 
consequently, we have estimated the intensity for 
that region. The true intensity could be quite differ- 
ent from what we have estimated, but the relative 
intensities of the lines in that region are probably 
good enough to recognize and assign the lines 
needed for calibration. That band is so weak that it 



is probably useless as a source of calibration for 
many workers, but it was included because it is 
based on frequency measurements. Since the lower 
energy level of the band is quite high, heating the 
absorption cell will make a large difference in the 
intensity. 

Please see Eqs. (3.24) to (3.29) for the definition 
of the various Herman-Wallis constants used in this 
work. In his intensity studies [5.218, 5.185] Toth has 
given Herman-Wallis terms for all the bands. For 
the most part we have used those constants in the 
intensity calculations for the tables and figures. 
Several exceptions were made, however. For the 
OS'O-Ol'O hot band we have used the same Herman- 
Wallis terms as given by Toth for the 02°0-00°0 
band. For the other transitions involving Av 2 =2, we 
have estimated the intensity constants to be used 
for the calculations because they were not included 
in Toth's tables. 

There seems to be no determination of a 
Herman-Wallis constant for the Av 2 = 1 transitions, 
therefore, we used a\ = 0.0016 which was estimated 
from the experimental spectrum. This value was es- 
timated from laboratory spectra that were com- 
pared with the calculated spectrum. The intensities 
of the lines in the bands in the 1835-1925 cm -1 re- 
gion were all calculated with £1 = - 0.0107 as mea- 
sured by Toth and Farmer [5.185]. The two 
Herman-Wallis constants measured by Boissy et al. 
[5.187] for the v x band, were used for all the Avi = 1 
transitions including hot bands and different isoto- 
pomers used in the tables. The 2vi + vi band also 
seemed to require a Herman-Wallis term although 
the spectra were better matched with a value of 
£2 = 0.3 x 10" 4 rather than the term given by Levy et 
al. [5.216]. No Herman-Wallis terms were used in 
the intensity calculations for the v\ — v$, v\ — 2vi, 2ps, 
and v\ + V2 bands and the hot bands that accompany 
them. 

3.4.3 CS 2 Only a few intensity measurements 
have been made on the vs band of CS 2 [5.375, 5.377, 
5.378, 5.379]. There seem to be no measurements of 
individually resolved rovibrational transitions, only 
integrated band intensities. For the purposes of this 
atlas we have calculated all hot band and isoto- 
pomer transitions with the same transition moment, 
0.27 debye, and with no Herman-Wallis terms. The 
integrated intensity for the 1500 cm -1 band region 
of CS 2 is 9.309 xlO -17 cm/molecule (2308 cm -2 
atm -1 at 296 K) as determined from adding all the 
line intensities including hot bands and different 
isotopes. This may be compared to the value 
9.55 x 10" 17 cm/molecule given by McKean et al. 
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p.j/oj, v.joxiu cm/moiecuie given Dy 

Robinson [5.375], 9.23 xlO" 17 cm/molecule given 
by Kiyama and Ozawa [5.377], and 9.13 x 10~ 17 cm/ 
molecule given by Person and Hall [5.379]. 

3.4.4 CO 

Fundamental Band 

Because of the large spacing between lines, CO 
was one of the first molecules for which individual 
line intensities were measured [5.255]. More recent 
measurements are tabulated in the review article 
by Smith et al. [5.2]. For the fundamental band 
near 2143 cm -1 it is difficult to determine which 
measurements are best but most of the recent mea- 
surements give integrated band intensities close to 
1.027xl0" 17 cm/molecule (276 cm -2 atm -1 at 
273.15 K). For this atlas we have used a dipole 
transition moment of 0.1073 D which gives a total 
band intensity of 1.03 xlO" 17 cm/molecule (276 
cm -2 atm -1 at 273.15 K or 255 cm -2 atm -1 at 296 
K). The changes in the matrix elements for the less 
abundant isotopic species were taken from the the- 
oretical calculations of Chackerian and Tipping 
[5.305]. The isotopic abundance was taken from 
Table 11. The uncertainty in the total band inten- 
sity seems to be about 3 percent. The intensity 
measurements of Chackerian et al. [5.300] 
(1.027 xlO -17 cm/molecule) are in agreement with 
other recent measurements for the main isotopic 
species and they also have studied the intensity of 
the weaker isotopic transitions. They gave transi- 
tion intensities calculated from the electric dipole 
moment function given by Chackerian and Tipping 
[5.305]. Their calculated intensities seem to include 
a weak Herman-Wallis effect but they do not give 
any explicit constants for easily calculating that ef- 
fect. A calculation that duplicates their tempera- 
ture and dipole derivative agrees with their values 
to within one percent (for / < 35) even without in- 
cluding a Herman-Wallis effect. 

In an earlier paper [5.282] Tipping gave calcu- 
lated values for the Herman-Wallis constants but 
they were small enough to ignore for this atlas. 
Bouanich [5.310] has also given calculated values 
for the Herman-Wallis constants. Apparently there 
are only two experimental determinations of 
Herman-Wallis constants that have been reported 
for the fundamental band of CO [5.261, 5.318]. The 
earlier work is subject to question because of the 
sensitivity to temperature errors. The most recent 
measurement was reported after this work was 



done and would change the intensities by no more 
than one percent. The intensity calculations used 
for the atlas did not include a Herman-Wallis 
effect for the fundamental band. 

First Overtone 

Fewer intensity measurements have been made 
on individual lines of the overtone band near 4260 
cm -1 . A good average of the more recent values 
seems to give an integrated band intensity of 
7.78 xlO -20 cm/molecule (2.09 cm -2 atm -1 at 
273.15 K) with an uncertainty of about 6 percent. 
This is equivalent to a transition dipole matrix ele- 
ment of 0.0066 D. 

The Herman-Wallis effect is significant for the 
first overtone band and has been included in the 
intensity calculation for the atlas. The constants 
given by Tipping [5.282], G = 0.005 and C 2 = 0.000 
034 [see Eq. (3.25)], were used for the calculation 
of the intensities of the first overtone transitions. 
Those calculated Herman-Wallis constants were in 
agreement with three experimental determinations 
[5.265, 5.271, 5.295] as well as with the calculation 
of Toth et al. [5.265]. 

3.4.5 NO A transition moment of 0.00412 D 
was used to calculate the line intensities given in 
the tables. This gives an integrated band intensity 
of about 5.04 x 10" 18 cm/molecule (125 cm" 2 atm" 1 
at 296 K) which agrees with the measurements of 
King and Crawford [5.330], Mandin et al. [5.359], 
and Holland et al. [5.368]. Some earlier measure- 
ments [5.326, 5.344, 5.347, 5.364, 5.367] and even 
some recent measurements [5.374] indicate that 
the intensity might be more like 4.44 to 4.64 x 10 ~ 18 
cm/molecule (110 to 115 cm" 2 atm" 1 at 296 K) so 
the intensities given in the tables may be too large 
by about 10 percent. We have taken the higher in- 
tensity value because NO is prone to having impu- 
rities that are hard to remove, thus giving low 
intensity readings. 

3.5 Other Line Parameters 

3.5.1 Lineshape and Pressure Broadening 

Smith et al. [5.2] have given a good discussion of 
lineshapes and the determination of pressure 
broadening coefficients. We shall give here only a 
brief description to ensure that the reader under- 
stands the principal equations describing these ef- 
fects. For a more complete understanding of the 
subject one should refer to the work of Smith et al. 
[5.2] and the papers to which they refer. 
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For a static gas in a field-free environment at 
pressures below one atmosphere, there are four 
major factors that might contribute to the shapes of 
infrared absorption lines, 1) lifetime broadening, 2) 
Doppler broadening, 3) pressure broadening, and 
4) collisional narrowing. 

Lifetime Broadening 
The lifetime broadening of a state n is given by 

y„=h/2irr„ , 

where y„ is the half-width of the state and t„ is the 
lifetime of that state. For a transition i-*j between 
two states, i and;, the transition linewidth will be 
given by the lifetime of both upper and lower 
states, 

m = (l/n + l/ij)*/2ir. 

For most stable molecules in the ground electronic 
state the radiative lifetime in a given rovibrational 
state will be on the order of 1 ms or more which 
gives a linewidth of about 1 X 10" 8 cm -1 (0.16 kHz) 
or less. This is much smaller than the Doppler 
width even at temperatures on the order of 5 K. 
For Doppler-free measurements the effect of life- 
time broadening could be important even for stable 
molecules although various instrumental effects, 
such as beam width or beam collimation, usually 
limit the effective linewidth. Transit time broaden- 
ing is a variation of lifetime broadening where the 
lifetime is the time that the molecule is in the light- 
beam. 

Doppler Broadening 

The Doppler width is the result of the random 
motion of the molecules in a gas sample and is 
given by 



7b = 3.581 xl(T 7 Kr/M) 



1/2 



(3.30) 



where v is the wavenumber or frequency of the 
transition, T is the temperature of the gas in kelvin, 
M is the relative molecular mass of the molecule in 
atomic mass units, and td is half the width of the 
transition at half the intensity (half-width at half 
height or HWHH). Since spectrometer resolution 
is often expressed in terms of the full width at half 
the line height (FWHH), the Doppler width is 
sometimes given by 2yo. 



For most gases at ambient temperatures and 
pressures below 2 kPa (15 Torr), the true lineshape 
(that is, the lineshape that would be observed with 
an instrument with infinite resolution) is domi- 
nated by the Doppler effect. For some molecules 
the effect of collisional narrowing, see below, is 
also important in that pressure regime. The 
Doppler effect gives a Gaussian lineshape which is 
described by the function 



/(^)={[(ln2)M" 2 / 7D } 
exp{-(ln2)[(^-M ) )/7D] 2 } 



(3.31) 



where u is the frequency (or wavenumber) of the 
center of the line. 

In Eq. (3.31) the Gaussian shape function f{v) 
has been normalized so that 



f /(„)d„ = l. (3.32) 



Consequently, since 

k(u) = Sf(u), 



(3.33) 



integrating both sides of Eq. (3.33) over all fre- 
quency space gives 

f k(v)dv= f Sf(v)dv = S. (3.34) 

J — 00 J — oo 



Here k{v) is the absorption coefficient at fre- 
quency v and 5 is the integrated line intensity given 
earlier in Eq. (3.9). 

One important characteristic of the Gaussian 
lineshape is the small wing absorption due to the 
exponential reduction in absorption as one gets far- 
ther from the center of the line. As pointed out by 
Korb et al. [5.263], accurate intensity measure- 
ments are more easily made when the wing absorp- 
tion is small. Another characteristic of the 
Gaussian shape is the bluntness at the center of the 
line. 

Pressure Broadening 

Pressure broadening gives rise to lines with a 
Lorentzian lineshape for which the normalized 
shape function has the form 

m = {ydir)l[{v-vtf + yll (3.35) 
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where yt is half the width of the line at the half 
intensity point. When pressure broadening is the 
dominant effect determining the lineshape, then 
Eq. (3.35) must be used in Eqs. (3.33) and (3.34). 

Pressure broadening is an additive effect, there- 
fore, the broadening of each gas in a mixture de- 
pends only on the partial pressure, P x , of that gas. 
The total broadening in a mixture is the sum of the 
broadening of each gas. Thus, for gases a and b 
with broadening coefficients c a and Cb, the total 
pressure broadened width will be yL=c a Pa+CbPb. 
The broadening coefficient is unique to each ab- 
sorbing molecule and to each collision partner. The 
broadening coefficient is a function of the tempera- 
ture as one might expect since the average collision 
velocity changes with temperature. The broadening 
coefficient is also different for each rotational tran- 
sition although the changes are systematic with the 
rotational quantum numbers. 

The pressure broadening coefficients generally 
have values in the range of 400 to 1300 MHz/Pa (3 
to 10 MHz/Torr) and are greatest for molecules 
that have large dipole moments. The pressure 
broadening coefficients are generally smallest for 
those rotational energy levels that are at relatively 
high energy. 

In contrast to Gaussian shaped lines, the 
Lorentzian shaped lines are sharper at the line 
center but have very extensive wings. For very 
strong lines it is possible to have significant absorp- 
tion intensity twenty half-widths from the line cen- 
ter. This is noticeable in the CO atlas where there 
is only a small amount of pressure broadening and 
yet the strong lines have very noticeable wings. If 
there were no pressure broadening, the wings of 
the strong lines would not be so prominent. 

Voigt Profile 

In many cases the lineshape is determined both 
by the Doppler effect and by the effects of pressure 
broadening. In such cases the lineshape is more ac- 
curately given by a Voigt profile which is a convolu- 
tion of the Gaussian and Lorentzian profiles. 
There is no good single closed-form expression for 
the Voigt profile. Rather, the Voigt shape function 
is given by the integral expression 

f(v) = (Bx/v) f [exp(-y 2 )V[x 2 + (z- y y]<iy 

J — to 

(3.36) 

and various approximations to that integral. Equa- 
tion (3.36) has been simplified by using 



and 



B = (l/- m )[(ln2)/7r] 1/2 , 
* = (rL/rD)(ln2)" 2 , 

z = [(*-n,)/-tt>](ln2) ,/2 . 



There are a number of good computer programs 
[3.17, 3.18] for evaluating the Voigt shape function. 

Collisional Narrowing 

Collisional narrowing or Dicke narrowing [3.19- 
3.21] has only recently been measured for a few 
molecules, but its effects have been observed for 
NO [5.370]. Collisional narrowing has the effect of 
reducing the size of the Gaussian linewidth yo. In 
other words, it makes the line appear to have a 
Doppler width that is smaller than that calculated 
by Eq. (3.30). For NO, the effective Gaussian 
width is about 9 percent smaller than the Doppler 
width at a pressure of 6.6 kPa (50 Torr). Since col- 
lisional narrowing is primarily a kinetic collisional 
effect (sometimes described in terms of hard and 
soft collision models), it is not expected to depend 
on the rotational or vibrational levels involved in 
the transition. There is, however, a weak depen- 
dence on the transition assignment as shown by 
Pine and Looney's work on HC1 and HF [3.22]. 

3.5.2 Estimating Peak Intensities 

From tables of line intensities, such as are given 
in this atlas, it is possible to estimate the pressure- 
pathlength product needed to obtain a spectrum 
with adequate intensity. For this purpose we con- 
sider the two limiting cases of Doppler (or Gaus- 
sian) shaped lines at low pressures and the 
pressure broadened Lorentzian lines for high pres- 
sures. 

If the effect of pressure broadening is negligible 
and the lineshape is determined by the Doppler 
width rather than the spectrometer slit function, 
the percent transmission at the center of a line is 
given by 

% transmission = 100 exp( - CSlp /y ), (3.37) 

where C is 1.1494 x 10 14 if 5 is the intensity given in 
the tables (in units of cm/molecule), y is the 
Doppler width [given by Eq. (3.30)], / is the 
pathlength in centimeters, and p is the pressure in 
pascals (Pa). (If p is measured in Torr, then 
C = 1.5324 x 10 16 .) Note that C has been evaluated 
for a temperature of 296 K, which is the 
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temperature for which the intensities have been 
calculated in this atlas. 

If the lineshape is dominated by pressure broad- 
ening, the percent transmission is again given by 
Eq. (3.37), but y should be the pressure-broadened 
linewidth, yt, and C will be 7.789 x 10 13 if p is in Pa 
(1.038xl0 15 if/?isinTorr). 

For the same linewidth and integrated intensity a 
Doppler-broadened line will have a peak intensity 
1.476 times greater than a pressure-broadened line, 
even though the pressure broadened line is 
sharper. 

For intermediate pressures the peak intensity 
can be more accurately estimated by modifying Eq. 
(3.37) so as to add the approximate contribution of 
both shapes according to 

% transmission = 100 exp{- S/p[(1.1494 X 10 14 td 

+ 7.789 xl0 13 7t)/(->£ + yD]} (3.38) 

Equation (3.38) reduces to Eq. (3.37) when either 
yo or 7t dominates the lineshape. 

The peak intensity observed with an instrument 
that introduces any instrumental broadening will, 
of course, be smaller than that calculated with ei- 
ther Eq. (3.37) or (3.38). As a rule of thumb, the 
peak intensity will be diminished by more than the 
ratio •y/'ys, where y is the true linewidth and ys is 
the width of the instrumental resolution function. 
Thus, for an instrument with a resolution function 
that is ten times greater than the true linewidth, an 
absorption line will appear at least ten times 
weaker than what is calculated by either Eq. (3.37) 
or Eq. (3.38), provided the line is not saturated. 

3.5.3 Pressure Induced Lineshifts A good fre- 
quency calibration standard is one whose frequency 
is not changed as the measurement conditions are 
varied. Some of the absorption lines given in this 
atlas are weak enough to require that either long 
pathlengths or moderate pressures, 0.3 to 1.3 Pa (2 
to 10 Torr), be used. It is important that one recog- 
nize the additional calibration uncertainty intro- 
duced by using pressures that are too high. 
Although pressure induced shifts in the frequency 
of the absorption lines are poorly understood and 
measurements are few and not very accurate, it is 
possible to estimate the approximate effect of mod- 
erate pressures on the frequencies of these calibra- 
tion standards. 

As a general rule, pressure induced frequency 
shifts are at least an order of magnitude smaller 
than pressure induced broadening. For the over- 
tone of CO, Pollock et al. [5.304] found a pressure 



shift (due to CO) on the order of -2±1.5 kHz/Pa 
(-0.3 ±0.2 MHz/Torr). Bouanich [5.302] found 
pressure shifts on the order of -1.2 ±0.3 kHz/Pa 
(-0.1 6 ±0.04 MHz/Torr). He seems to have ob- 
served a significant rotational dependence, but we 
only quote an average value. For the overtone of 
NO, Pine et al. [5.370] found a self-induced pres- 
sure shift of -1.1 ±0.3 kHz/Pa (-0.15 ±0.04 
MHz/Torr). They found that there may be a weak 
rotational dependence but it was obscured by ex- 
perimental error. 

For N 2 an attempt to measure the self-induced 
pressure shift in the rotational spectrum resulted in 
an estimate of an upper limit of 0.75 kHz/Pa 
[5.201]. The pressure shift of several lines of N 2 
were measured near 4500 cm -1 and an average 
value of -1.2 ±0.3 kHz/Pa (-0.16 ±0.04 MHz/ 
Torr) was found [5.219]. The pressure shift caused 
by N 2 and 2 on the i> 3 band of N 2 near 1280 cm -1 
was measured by Varanasi and Chudamani [5.238]. 
Their diode laser measurements gave an average 
value of about -0.6 ±0.15 kHz/Pa (-0.08 MHz/ 
Torr) for several lines between 1250 and 1300 
cm -1 . 

More pressure shift measurements seem to have 
been made on OCS than on any of the other 
molecules in this atlas. The earlier measurements 
indicated that the pressure shift was on the order 
of 4.0 kHz/Pa (0.5 MHz/Torr) or less [5.94]. Later 
this estimate was improved by new measurements 
that gave an average value of —0.4 ±1.5 kHz/Pa 
(-0.05 ±0.20 MHz/Torr) [5.120]. The most recent 
measurements at about 1000 cm -1 [5.133] indicate 
that the self-shift is -0.37 ±0.04 kHz/Pa 
(-0.049 ±0.005 MHz/Torr). Recent measurements 
of the self-shift of microwave transitions [5.119] 
found that the shift was too small to measure, 
0.000 ±0.04 kHz/Pa or 0.000 ±0.006 MHz/Torr. 

Kou and Guelachvili [3.23] have recently mea- 
sured the self-induced pressure shift for the C0 2 
laser lines near 1000 cm -1 . They found the shift to 
be on the order of -1.05 ±0.2 kHz/Pa 
(-0.14 ±0.03 MHz/Torr) with no evidence of a/- 
dependence. One might expect that the self- 
induced pressure-shift of C0 2 would be similar to 
that of the molecules used in the present compila- 
tion. 

Since the pressure shift is a shift in the energy 
levels, it is likely to be greater as one goes to higher 
energy levels. As a first approximation one can 
probably assume that the shift is proportional to 
the frequency. It is also possible that for poly- 
atomic molecules, the pressure shift may depend, 
to a large extent, on the vibrational mode involved. 
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For purposes of estimating the maximum error 
that could be introduced in a calibration measure- 
ment using any of the data given in this atlas, one 
should treat the frequencies given in this work as 
applying for a pressure below 130 Pa (1 Torr). For 
each increase in pressure of 130 Pa (1 Torr) the 
uncertainty in the frequency for N2O, NO, CO, and 
CS2 should be increased by 0.3 MHz unless the 
pressure shift is added to the frequencies given in 
these tables or unless more accurate values for the 
pressure shift become available. For OCS the un- 
certainty should be increased by about 0.1 MHz/ 
Torr. Note that the pressure shift seems to be 
negative for all the molecules in this atlas. 
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4. Error Analysis 

The standards presented in this book are based 
on data derived in large part from frequency mea- 
surements because such measurements are less 
prone to systematic errors. In assessing the uncer- 
tainties that should be assigned to the recom- 
mended calibration frequencies, we must consider 
five major factors that may contribute to error in 
the determination of the line frequencies. 

1. The accuracy of the calibration source. 

2. The uncertainty in locating the center of the ab- 
sorption line. 

3. The accuracy of transferring the calibration to 
the line center measurement. 

4. Errors caused by environmental effects. 

5. Model errors, or uncertainties in the applica- 
tion of least-squares techniques to obtain the 
best estimate of the correct line frequency. 

In succeeding paragraphs of this section each of 
these factors will be examined in relation to the cal- 
ibration frequencies recommended in this book. 



4.1 Accuracy of the Calibration Source 

Much of the primary infrared data used for these 
tables originated at the NIST laboratory in Boulder, 
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Colorado. Two other laboratories have reported 
heterodyne frequency measurements on the 2in 
band of OCS, [5.88, 5.126]. The most important of 
these were two saturated absorption measurements 
made by Fayt et al. [5.126] which considerably re- 
duce the uncertainty in the frequencies for 2vi. The 
measurement of the laser transitions of N 2 by 
Whitford et al. [5.183] gave accurate frequencies 
for the separation of the v\ and vi states which 
were useful for the tables between 880 and 980 
cm -1 . Recently frequency measurements have been 
made by Urban and coworkers [5.321, 5.137a]. 

The NIST measurements used well character- 
ized CO2 lasers that were virtually identical to the 
lasers used by Petersen et al. [2.9, 2.11, 2.12] in the 
last major determination of the CO2 laser frequen- 
cies. Since much of the laser technology that went 
into the determination of the CO2 laser frequencies 
was developed in the same laboratory where the 
heterodyne measurements were made, one can be 
sure of the accuracy of the calibration source. 
These lasers have been described in Sec. 2. The 
CO2 laser frequencies are good to at least ±0.05 
MHz and are not a significant source of error for 
the Doppler limited heterodyne measurements. 

4.2 Uncertainty in Locating the Line Center 

Two major factors contribute to errors in the lo- 
cation or determination of absorption line centers: 
the signal-to-noise ratio, and the slope or other ir- 
regularities in the background. 

For the infrared measurements the slope or ir- 
regular background in the radiation being absorbed 
was only a minor contributor to error since its ef- 
fect could be practically eliminated, as described in 
Sec. 2. 

In some cases imperfectly resolved or overlap- 
ping lines were measured. In cases where the lines 
are close doublets, the measurement represents the 
center of gravity of the doublets, and the least- 
squares analysis took that into account. If the lines 
were partially resolved, that was taken into account 
in assigning the uncertainty for the measurement. 
In most cases partially resolved lines that have un- 
equal intensities should not be used unless a very 
generous allowance is made for the uncertainty 
since the resulting error is not random, but is in a 
particular direction. OCS, CO, and CS2 are good 
molecules to use for standards because they have 
no quadrupole fine structure. N2O is not as good 
because all of its transitions have a small quadru- 
pole splitting. Since the splitting for N2O is never 
greater than 4 MHz and diminishes rapidly with 
increasing /-values and since the Doppler full- 



width of N 2 at 500 cm" 1 is 28 MHz, the splitting 
should not affect the use of N2O for most infrared 
calibration purposes. 

For each infrared heterodyne frequency mea- 
surement two line centers must be measured: the 
absorption line center, and the center of the differ- 
ence frequency (between the TDL and the local 
oscillator). The location of the absorption line cen- 
ter is determined by the technique used to lock the 
TDL to the absorption line as described in Sec. 2. 
The uncertainty assigned to the line-center lock is 
given by half the Doppler-width of the line divided 
by the signal-to-noise ratio of the derivative signal 
for the line. Because of modulation broadening the 
linewidth is slightly greater than the Doppler-width 
but this approximation to the lock uncertainty is 
adequate. The lock error is random and is reflected 
in the statistical analysis of the least-squares fit of 
many measurements. 

More important than the uncertainty in locking 
to the line center is the uncertainty in the difference 
frequency measurements. That uncertainty is given 
by one-tenth the beat note linewidth, or by half the 
beat note width divided by the signal-to-noise ratio, 
whichever is larger. We believe this covers both ran- 
dom and residual systematic errors for a single mea- 
surement. In most of our measurements, this has 
been the predominant uncertainty and the lock un- 
certainty has been negligible by comparison. 

One source of error in the heterodyne measure- 
ments that has no direct counterpart in wavelength 
measurements is due to frequency-dependent dif- 
ferences in the transmission or amplification com- 
ponents involved in the heterodyne frequency 
measurements. This frequency dependence some- 
times presented an additional distortion to the het- 
erodyne frequency lineshape as displayed by the 
spectrum analyzer. Fortunately this error, which 
would be systematic for all measurements of the 
same heterodyne frequency, becomes randomized if 
enough different lines are measured. In most cases 
this problem was recognized and was taken into ac- 
count in estimating the uncertainty of each mea- 
surement. 

Examination of our results over the past decade 
indicates that the procedure outlined above for as- 
signing the uncertainty produced values which 
turned out to be close to a 2 or uncertainty. 

For the microwave and sub-millimeter wave 
measurements used to prepare these tables, the 
uncertainties given in the literature are generally 
accurate enough although we have increased the 
uncertainties slightly in a few cases. There may be 
some systematic error due to incomplete 
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modulation when Stark modulation was used, but 
this is compensated by the much smaller 
linewidths, and the smaller frequency dependent 
irregularities, compared to the infrared measure- 
ments. 

43 Accuracy of Transferring the Calibration to 
the Line Center Measurement 

The primary difference in the reliability of wave- 
length and frequency measurements is in the accu- 
racy of transferring the calibration to the mea- 
surement. In frequency measurements the accuracy 
of the measurements is not affected by beam di- 
mensions, by wavefront mismatch or other mis- 
alignments, or by refractive index effects. The 
frequency of the radiation is always the same no 
matter what the medium may be, or how it is mea- 
sured. Modern electronics excel at counting, and 
that is how frequency is determined. With fre- 
quency measurements, if a measurement can be 
made at all, the uncertainty must come from the 
four other sources of error discussed in this chap- 
ter. 

On the other hand, wavelength comparisons are 
susceptible to many different wavelength-depen- 
dent errors, errors due to misalignment, errors due 
to differences in the ratio of the dimensions of the 
optical elements (or of the optical beam), to the 
wavelength. The best wavelength measurements 
use a calibration that is nearly the same wavelength 
and intensity as the feature to be calibrated. Each 
type of wavelength measurement has its own pecu- 
liarities. 

The present tables are based in part on higher 
order centrifugal distortion constants determined 
to some extent by FTS measurements which are 
essentially wavelength measurements. These mea- 
surements were all internally calibrated by means 
of lines whose frequency could be determined by 
heterodyne frequency measurements. These cali- 
bration features were always within the approxi- 
mately 120 cm -1 band-pass of the FTS mea- 
surements. For the FTS measurements there may 
be a phase error which is different for weak fea- 
tures and for strong features and that is one reason 
why greater emphasis was placed on the use of fre- 
quency measurements of pure rotational transi- 
tions to determine the rotational constants 
wherever possible. In no case were FTS measure- 
ments used to determine the band centers or vibra- 
tional levels for the lines recommended as 
frequency standards. 



4.4 Errors Caused by Environmental Effects 

For closed shell molecules the only significant 
environmental effects are due to pressure shifts 
and electric field effects such as either the ac or dc 
Stark effects. As long as the radiation field (for 
absorption spectra) is too low to give saturated ab- 
sorption effects, the ac Stark effect can be ignored 
in its effect on the absorption line center. Even the 
dc Stark effect can be ignored for most work since 
stray electric fields are generally too small to give 
noticeable Stark shifts. 

We believe that even in the case of the saturated 
absorption measurements on the 2-0 band of CO, 
the shift due to the ac Stark effect will be smaller 
than a tenth of the linewidth, or less than 0.2 MHz. 

The ac Stark effect may have a small effect on 
the frequency of the CO2 lasers but, since those 
lasers have the same characteristics as the lasers 
used in the original measurements against the ce- 
sium frequency standard, the effect of the shift will 
already be included in the frequency assigned to 
the laser lines. In the case of the CO transfer oscil- 
lators, the CO frequencies are measured at the 
same time as the beat note frequency measurement 
and any ac shift is included in the measurement. 

On the other hand, the effect of pressure in- 
duced shifts is potentially significant and deserves 
serious consideration. Pressure-induced shifts in 
infrared spectra have not been extensively studied, 
so there is no experimentally confirmed theory that 
one can use to calculate the pressure shifts to be 
expected for much of the data given in the present 
tables. All of the pressure shift measurements 
seem to indicate that the shift for the transitions 
given in these tables may be on the order of 
-2.2kHz/Pa (-0.3 MHz/Torr) or less. Only in the 
case of the measurements given by Vanek et al. 
[9.133] is the shift of -0.37 ±0.04 kHz/Pa (-49 ±5 
kHz/Torr) reliably given for the 2vz band of OCS. 
Since neither the frequency dependence nor the 
rotational or vibrational dependence of the pres- 
sure shift is known, this remains one of the most 
important uncertainties in the application of these 
tables to real measurements. 

4.5 Uncertainties in the Application of Least- 
Squares Techniques 

Because individual measurements of infrared ab- 
sorption lines by heterodyne measurement tech- 
niques are not very precise (uncertainties of the 
order of 5 to 10 MHz for some regions) and only a 
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small number of transitions can be measured, it is 
necessary to use least-squares techniques to com- 
bine all available measurements to yield calculated 
transitions that are considerably more accurate 
than any one measurement. Such fitting techniques 
can be no more accurate than the equations used 
to describe the transitions (the Hamiltonian). The 
molecules and the particular transitions in these 
tables were chosen in part because of the reliability 
of the Hamiltonian as shown by extensive studies 
reported in the literature. Although certain inter- 
actions, such as Fermi resonance, affect the 
constants used in the effective Hamiltonian, such 
resonances can be ignored, provided that measure- 
ments are available for a wide range of rotational 
levels for each vibrational state. In the analysis 
used in this book the effect of /-type resonance was 
included because it represents the largest 
resonance effect and it can be reliably estimated, 
thereby giving the correct functional form to the 
centrifugal distortion constants. The details of the 
fitting procedure and Hamiltonian are given in 
Sec. 3. 

The variance-covariance matrix determined by 
the least-squares fit gives a reliable estimate of the 
uncertainties of the transitions in the range of rota- 
tional levels for which there are good measure- 
ments. The calculated uncertainties get large quite 
rapidly for transitions extrapolated beyond the 
range of measured energy levels, but the reliability 
of those calculated uncertainties deteriorates even 
more quickly. For that reason we have terminated 
the recommended calibration standard indication 
at the highest /-value for which there are good 
measurements. Higher transitions are given in the 
tables but their accuracy is less certain. 

4.6 Summary 

Of the five sources of error identified above, the 
first (the accuracy of the calibration source) does 
not contribute significantly to the uncertainties in 
the calibration frequencies given in these tables. 
For frequency measurements the second and third 
source (uncertainty in locating the line center and 
accuracy of transferring the calibration) will appear 
as random errors and so will be given by the statis- 
tical analysis of the least-squares fit of the mea- 
surements. For the molecules and bands 
represented in these tables the uncertainties con- 
tributed by model errors are small and likely to 
show as deviations that are included in the statisti- 
cal analysis. 



In conclusion we think that the uncertainties 
given by the statistical analysis are adequate to de- 
scribe the errors in the frequency measurements 
that might arise from all causes except errors due 
to pressure-induced frequency shifts. Most of the 
heterodyne measurements were made at low pres- 
sures but some measurements of weaker transitions 
were made at pressures as great as 1200 Pa (9 
Torr). To allow an extra margin of error due to 
pressure shifts, the uncertainties assigned to the 
heterodyne frequency measurement data used in 
the fit were about twice as large as the rms devia- 
tion. Primarily, this had the effect of increasing the 
uncertainty in the vibrational energy levels. 

5. Bibliography 

We attempt here to list all references to papers 
giving infrared or microwave frequency (or wave- 
number) measurements, as well as lineshape and 
intensity measurements that are relevant to this at- 
las. Papers involving bands not included in this atlas 
may be missing from this bibliography. Some papers 
are not included if they involve foreign gas broaden- 
ing measurements only. For a more complete listing 
of pressure broadening papers see the review given 
by Smith et al. [5.2]. When completely superseded, 
some of the older papers may not appear in this bib- 
liography but they can be found referenced in the 
more recent papers. The references are grouped by 
molecule and arranged in chronological order for 
each molecule. At the end of each reference is a list 
of initials that indicate the subject matter covered 
by the reference as follows: 

F —frequency measurements, 

FB— foreign gas broadening measurements, 

I —intensity measurements, 
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6. Atlas and Wavenumber Tables 

6.1 Scope of the Atlas 

The range of wavenumber coverage of the 
molecular bands is indicated below for the five 
molecules selected for use. 

CO Atlas -1948 to 2275 cm -1 
and 4071 to 4352 cm -1 



OCS Atlas 
OCS Atlas 
OCS Atlas 
OCS Atlas 
OCS Atlas 
OCS Atlas 
OCS Atlas 
OCS Atlas 
OCS Atlas 
OCS Atlas 
N 2 Atlas - 
N 2 Atlas- 
N 2 Atlas - 
N 2 Atlas - 
N 2 Atlas - 
N 2 Atlas - 
N 2 Atlas - 



-486 to 567 cm" 1 
-812 to 890 cm" 1 
-1000 to 1095 cm" 1 
-1650 to 1739 cm" 1 
-1832 to 1934 cm -1 
-1970 to 2141 cm" 1 
-2510 to 2600 cm -1 
-2693 to 2763 cm" 1 
-2862 to 2970 cm -1 
-3065 to 3120 cm" 1 
-523 to 659 cm" 1 
-880 to 1087 cm -1 
-1105 to 1345 cm" 1 
-1820 to 1925 cm" 1 
-2140 to 2269 cm" 1 
-2400 to 2607 cm -1 
-2725 to 2842 cm -1 



CS 2 Atlas -1460 to 1551 cm -1 



NO Atlas -1741 to 1940 cm" 1 

6.1.1 Description Throughout the wavenum- 
ber tables the lines that are suitable for use as 
wavenumber standards are indicated by an asterisk 
( * ) following the wavenumber and its uncertainty. 
In assigning the asterisks no consideration was 
given to problems related to overlapping with other 
transitions. The tables list nearby lines that may 
cause problems with overlapping. The user must 
exercise judgment in determining if such overlap- 
ping will impair the accuracy of the measurement. 
The asterisk only certifies the accuracy of the line 
position in the hypothetical absence of any other 
nearby lines. Obviously, the resolution of the in- 
strumentation being used will determine if a 
nearby line might invalidate the accuracy of a cali- 
bration line. 

The uncertainties of the wavenumbers are given 
in parentheses after the wavenumbers in those 
cases where there is reason to believe that a good 
estimate of the uncertainty can be made. Even so, 
the uncertainty in the lines not designated as cali- 
bration lines should be taken with some degree of 
skepticism. 



The uncertainties given in the tables are twice 
the estimated standard error as calculated from the 
variance-covariance matrix given by the least- 
squares fit that determined the constants used to 
calculate the wavenumbers. The uncertainties refer 
to the accuracy of each individual transition. In 
general, the wavenumber separation of two nearby 
lines for the same vibrational transition of the same 
molecular species will be given more accurately 
than the uncertainty given in parentheses might 
lead one to believe. That is because the relative 
differences between the rotational energy levels 
are usually known more accurately than the differ- 
ences in the vibrational energy levels. 

On the other hand, the separation of two lines 
that are due to absorption from two different iso- 
topic species is probably known no more accurately 
than the uncertainty (given in parentheses) would 
lead us to believe. 

The wavenumbers given in the tables are calcu- 
lated wavenumbers because they are more reliable 
than individual line measurements and the uncer- 
tainties in the calculated wavenumbers can be ac- 
curately estimated. 

The wavenumber tables also contain a column 
for the intensity estimated for each transition at a 
temperature of 296 K. The format for the intensity 
values is the standard computer format consisting 
of a decimal value followed by the exponent (the 
power of ten multiplying the decimal value). The 
intensities given in the wavenumber tables are rep- 
resented by S in Eq. (3.10) and are integrated line 
intensities rather than peak intensities. 

In the sections giving discussions and equations 
on intensity calculations and on pressure broaden- 
ing (Sec. 3), equations are given for estimating the 
appearance of the spectrum for different experi- 
mental conditions. In particular, Eq. (3.37) can be 
used to estimate the percent transmission at the 
center of a line under different conditions of pres- 
sure and pathlength. As an aid in calculating inten- 
sities at different temperatures the tables for OCS 
and N 2 also contain a column giving the separa- 
tion (in cm -1 ) of the lower state energy level from 
the ground state. The units given in the tables can 
be converted to the more common units of cm -2 
aim -1 at 296 K by multiplying by 2.479 xlO 19 . To 
convert to intensities at some other temperature 
one should refer to Eq. (3.10). Smith et al. [5.2] give 
a table for converting to other units. 

The intensity values given in this work are only 
given as an aid in estimating the appearance of the 
spectrum, they should not be treated as well deter- 
mined values. The intensities given for weak lines 
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and especially for the rarer isotopes may be in er- 
ror by 50 percent or more. 

The spectral illustrations were actually calcu- 
lated spectra rather than reproductions of real 
measurements. This gave us more flexibility in 
choosing effective pressures and pathlengths that 
seemed most appropriate to illustrate even the 
weak lines. As with any digitized spectrum, regard- 
less of whether it is calculated or measured, the 
peak intensity of sharp lines may show some irregu- 
larity depending on whether the true peak falls on 
a digitized point or slightly misses it. The spectra 
were plotted with a digitizing interval of about 
0.0005 to 0.001 cm -1 . Close doublets that should 
have the same intensity may show slight intensity 
differences because of this digitizing effect. 

Comparison with real spectra measured in our 
own laboratory or illustrated in published works 
(such as Refs. [5.101, 5.102, 5.307]) showed that 
the spectra given in this work are adequate for 
identifying the calibration lines. Some weak transi- 
tions may be absent from the calculated spectrum 
even though they might be found in a real spec- 
trum of comparable pressure and pathlength. Cer- 
tainly, absorption due to common impurities such 
as H2O or CO2 will not be found in these spectra. 

The spectra used in the illustrations for CO were 
calculated for infinite resolution but the lines were 
given widths dictated by the Doppler width of the 
line convolved with the pressure broadened width. 
For the CO spectra the shape of the weak lines is 
dominated by the Doppler width which is 0.0047 
cm -1 (FWHH) at 2000 cm -1 and 0.0093 cm" 1 at 
4000 cm -1 . For the conditions chosen for the atlas 
illustrations, the strong CO lines are much broader 
than the weak ones and show pronounced shoul- 
ders due to the effect of even a veiy small pressure 
broadening. 

The figures used for illustrating the OCS and 
N2O spectra were calculated for spectrometer reso- 
lutions on the order of 0.003 cm -1 . Again the ex- 
pected Doppler and pressure-broadened 
lineshapes were used in calculating the spectra. 
The atlas for OCS and N 2 is divided into sections 
according to the vibrational transitions involved. At 
the beginning of each section the parameters (slit 
width, dipole derivative, Herman-Wallis constants) 
used in calculating the spectra are given as well as 
a key to the abbreviations used for the vibrational 
transitions in the atlas. 

In calculating the pressure-broadened width, a 
single value for the pressure broadening was used 
for all lines in a spectrum. It is well known that 



such an assumption is incorrect and therein lies 
one reason for the illustration to depart slightly 
from a true spectrum. 

The spectra overlap slightly in order to show the 
relationship of the lines near the ends of each 
panel. At the top of each spectrum fiducial marks 
indicate which lines are given in the accompanying 
wavenumber tables. Every fifth line is indicated 
with a darker and longer mark. In the tables, every 
fifth line is set apart by a following blank line. Even 
though the panels overlap, each line is identified by 
a fiducial mark only once. The next section gives a 
key to symbols in the atlas, followed by a sample 
spectral map with facing table. 
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6.2 Key to Symbols and Sample Spectral Map 
with Facing Table 



Band 


Isotopomer 


Vibrational 








transition 


A 


» 6 12 C? 2 S 


02 0-00°0 


X 






02 2e 0-00°0 


B 






(B'O-Ol 1 ^ 


C 






(tt'O-Ol^O 


D 






02 1-00°1 


E 






04 2 0-02 2e 


F 






04 2 0-02 2f 


G 






04 0-02°0 


H 






03'l-01 le l 


I 






03'l-01 lf l 


J 






05 3 0-03 3e 


K 






05 3 0-03 3f 


L 






05'0-03 lc 


M 






os'o-os^o 


P 


16QI 


2 C 34 S 


02°0-00°0 


Q 






(B'O-Ol^O 


R 






OS'O-Ol^O 


T 


IfiQl 


3 C 32 S 


02 0-00°0 



If the e or / designation is not specified (for 
/*0), then the transition is to either level, depend- 
ing on the selection rules and the change in the 
rotational quantum number, /. Spectra are given 
for a slitwidth of 0.002 cm -1 and a temperature of 
296 K. For the A, X, D, P, and T bands a transition 
moment of 0.0333 D was used and no Herman- 
Wallis constant was included in the intensity calcu- 
lation. For the other bands a transition moment of 
0.032 D was used with a Herman-Wallis constant 
of C, = 0.0019. 
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LINE WAVENUMBER(unc) LOWER STATE INTENSITY ASSIGNMENT LINE WAVENUMBER(unc) LOWER STATE INTENSITY ASSIGNMENT 






# 


(cm') 


(cm 1 ) 


(cm/molecule) 


ROTATION BAND 


# 


(cm" 1 ) 


(cm" 1 ) 


(cm/molecul* 


1 


1060.056 36( 6)* 


575.72 


0.739E-21 


R(16) 


C 


36 


1062.885 66(24) 


340.66 


0.130E-21 


2 


1060.176 77(10)* 


1574.38 


0.104E-22 


R(59) 


D 


37 


1062.935 247(14)* 


270.13 


0.373E-20 


3 


1060.214 07(24) 


249.29 


0.173E-21 


R(35) 


P 


38 


1062.963 36( 6)* 


632.52 


0.805E-21 


4 


1060.239490(13)* 


188.62 


0.464E-20 


R(30) 


A 


39 


1063.093 94(6)* 


632.64 


0.805E-21 


5 


1060.410 14(6)* 


582.57 


0.759E-21 


R(17) 


B 


40 


1063.165 09(9)* 


1078.78 


0.980E-22 


6 


1060.48640(6)* 


582.63 


0.758E-21 


R(17) 


C 


41 


1063.267 42(15) 


692.18 


0.329E-22 


7 


1060.656 05(24) 


263.53 


0.166E-21 


R(36) 


P 


42 


1063.335 56(24) 


357.27 


0.122E-21 


8 


1060.661 66( 9)* 


1055.59 


0.583E-22 


R(6) 


G 


43 


1063.389 156(14)* 


285.13 


0.356E-20 


9 


1060.685 510(13)* 


201.19 


0.451E-20 


R(31) 


A 


44 


1063.392 04(6)* 


642.27 


0.802E-21 


10 


1060.833 43( 6)* 


589.88 


0.775E-21 


R(18) 


B 


45 


1063.533 07( 6)* 


642.40 


0.801E-21 


11 


1060.917 75( 6)* 


589.95 


0.775E-21 


R(18) 


C 


46 


1063.585 80( 9)* 


1084.07 


0.103E-21 


12 


1061.076 43(9)* 


1058.44 


0.659E-22 


R(7) 


G 


47 


1063.703 96(15) 


704.08 


0.322E-22 


13 


1061.132 833(13)* 


214.17 


0.437E-20 


R(32) 


A 


48 


1063.786 81(24) 


374.27 


0.115E-21 


14 


1061.257 62(6)* 


597.60 


0.788E-21 


R(19) 


B 


49 


1063.82162(6)* 


652.42 


0.796E-21 


15 


1061.350 39(6)* 


597.68 


0.787E-21 


R(19) 


C 


50 


1063.844 398(14)* 


300.54 


0.339E-20 


16 


1061.388 43(16) 


648.43 


0.349E-22 


R(25) 


Q 


51 


1063.973 50(6)* 


652.56 


0.795E-21 


17 


1061.492 19(9)* 


1061.69 


0.731E-22 


R(8) 


G 


52 


1064.007 52( 9)* 


1089.77 


0.108E-21 


18 


1061.543 93(24) 


293.19 


0.151E-21 


R(38) 


P 


53 


1064.252 10( 6)* 


662.98 


0.787E-21 


19 


1061.581464(14)* 


227.55 


0.422E-20 


R(33) 


A 


54 


1064.300 978(14)* 


316.35 


0.322E-20 


20 


1061.682 71(6)* 


605.72 


0.797E-21 


R(20) 


B 


55 


1064.415 24(6)* 


663.13 


0.786E-21 


21 


1061.784 33(6)* 


605.81 


0.797E-21 


R(20) 


C 


56 


1064.430 25( 9)* 


1095.88 


0.112E-21 


22 


1061.90892(9)* 


1065.36 


0.800E-22 


R(9) 


G 


57 


1064.683 48( 6)* 


673.94 


0.775E-21 


23 


1061.989 85(24) 


308.62 


0.144E-21 


R(39) 


P 


58 


1064.693 35(24) 


409.47 


0.102E-21 


24 


1062.031 406(14)* 


241.34 


0.406E-20 


R(34) 


A 


59 


1064.758 903(14)* 


332.57 


0.305E-20 


25 


1062.10869(6)* 


614.25 


0.803E-21 


R(21) 


B 


60 


1064.854 00( 9)* 


1102.39 


0.115E-21 


26 


1062.219 56(6)* 


614.35 


0.803E-21 


R(21) 


C 


61 


1064.858 28( 6)* 


674.10 


0.775E-21 


27 


1062.326 65( 9)* 


1069.42 


0.864E-22 


R(10) 


G 


62 


1065.115 76(6)* 


685.31 


0.761E-21 


28 


1062.398 22(16) 


669.58 


0.341 E-22 


R(27) 


R 


63 


1065.148 66(24) 


427.66 


0.949E-22 


29 


1062.437 09(24) 


324.44 


0.137E-21 


R(40) 


P 


64 


1065.218 178(14)* 


349. 19 


0.289E-20 


30 


1062.482 665(14)* 


255.53 


0.390E-20 


R(35) 


A 


65 


1065.278 77( 9)* 


1109.30 


0.118E-21 


31 


1062.535 58( 6)* 


623.18 


0.806E-21 


R(22) 


B 


66 


1065.548 94( 6)* 


697.09 


0.746E-21 


32 


1062.652 89(15) 


680.52 


0.336E-22 


R(28) 


Q 


67 


1065.678 810(14)* 


366.22 


0.272E-20 


33 


1062.656 10( 6)* 


623.29 


0.805E-21 


R(22) 


c 


68 


1065.704 57(9)* 


1116.63 


0.121E-21 


34 


1062.745 37( 9)* 


1073.90 


0.925E-22 


R(ll) 


G 


69 


1065.748 25( 6)* 


697.27 


0.745E-21 


35 


1062.832 18(16) 


680.68 


0.336E-22 


R(28) 


R 


70 


1065.983 03( 6)* 


709.27 


0.728E-21 



R(41 

R(36 

R(23 

R(23 

R(12 

R(29 
R(42 
R(37 
R(24 
R(24 

R(13 
R(30 
R(43 
R(25 
R(38; 

R(25 
R(14 
R(26 
R(39 
R(26 

R(15 
R(27 
R(45 
R(40; 
R(16; 

R(27 
R(28 
R(46 
R(41 
R(17 

R(29 
R(42 
R(18 
R(29 
R(30; 



P 
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B 
C 
G 

R 
P 
A 
B 
C 

G 
R 
P 
B 

A 

C 
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B 
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C 

G 
B 
P 
A 
G 

C 
B 
P 
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G 

B 
A 
G 
C 
B 
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7. Appendix A: Useful Equations and 
Formulas 

Conversion factors 

dipole moment or transition moment: 
lD=3.335 64xlO- 30 Cm 

wavenumber: 

lcm" ' = 29 979.2458 MHz 

pressure: 
1 atm = 101 325. Pa (pascal) 
lTorr = 133.322 Pa 

absorption intensity: 
1 cm" 2 atm-1 at 296 K=4.033xl0- 20 cm/ 
molecule 

Useful equations 

Boltzmann factor: 

N t =Nj exp (-1.439 Eij/T), 

where Ni is the population of the *th level and Nj is 
the population of the/th level with energy differ- 
ence in cm -1 of Ejj at a temperature of T kelvins. 

Doppler half-width at half intensity, yo: 

yo = 3.581 X 10" 7 v (7/M ) 1/2 , 

where v is the frequency in the same units as jo, T 
is the temperature in kelvins, and M is the mass of 
the molecule in atomic mass units. 

Peak intensity of a Doppler shaped line: 

% transmission = 100 exp ( - 1.1494 X 10" Slplyo), 

where S is the intensity given in the tables, / is the 
pathlength in cm, and p is the pressure in Pa. 

Peak intensity of a Lorentzian line: 

% transmission = 100 exp ( - 0.7789 x 10 M Slp/yt), 

where S is the intensity given in the tables, jl is the 
Lorentzian half-width at half height, / is the path- 
length in cm, and p is the pressure in Pa. 
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8. Appendix B. Energy Level Diagrams for OCS and N 2 
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Information on NIST SP821 and Standard 
Reference Database 39 

NIST Special Publication 821 Wavenumber Cali- 
bration Tables from Heterodyne Frequency Mea- 
surements, by Arthur G. Maki and Joseph S. Wells, 
contains not only the spectral maps with facing 
tables, but also includes the text duplicated in this 
journal. A limited number of hard bound copies 
are still available from the authors, or NIST SP 821 
may be ordered from: 

Superintendent of Documents, 
U.S. Government Printing Office, 
Washington, DC 20402 

Some users are not interested in spectral maps 
and prefer to have only the tabular information, 
(for example to load files into a computer). For 
these users we offer an alternative of having only 
diskettes of the tables. We have included energy 
level diagrams of both OCS and N 2 in Appendix 
B for the convenience of these users. 

In preparing the tables for this atlas, far more 
transitions were calculated than could possibly be 
published. A complete list of transitions calculated 
for this work is available as NIST Standard Refer- 
ence Database 39, Wavelength Calibration Tables, 
in the form of four diskettes for a personal com- 
puter. These diskettes, written in ASCII, contain 
the following information: wavenumber (cm -1 ), un- 
certainty, lower state energy (cm -1 ), intensity (cm/ 
molecule), vibrational assignment, rotational 
assignment, date of data entry, isotopic species, 
and band designation according to the keys given in 
this book. The authors intend to update the data- 
base periodically. Persons interested in obtaining 
the most up-to-date calibration data should con- 
tact: 

Standard Reference Data, 

National Institute of Standards and Technology 

Bldg. 221/Room A320, 

Gaithersburg, Maryland 20899 

(301) 975-2208 

Fax (301) 926-0416 
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